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I.  INTRODUCTION 


"The  global  trade  'n  conventional  armaments  has  become  a  burning  political 
issue  in  recent  years,  owing  primarily  to  the  dramatic  increase  in  the  volume  of  the 
international  market  since  the  early  I'^'iVs."  [Ref.  1:  p.  1]  The  new  significance  in 
world  arms  sales  is  the  result  of  three  factors.  [Ref  2;  pp.  275-276]  The  first  is  the 
sheer  '.olume  of  arms  being  traded  and  the  increasing  quality  of  these  weapons.  It 
used  to  be  the  norm  that  countries  would  sell  only  obsolete,  unneeded  or  second  rate 
equipment.  Now  nations  are  exporting  state  of  the  art  weapon  systems  with 
unparalleled  accuracy  and  destructive  capacity.  The  second  factor  is  the  decline  of 
effectiveness  in  alternative  methods  for  the  superpowers  to  influence  other  countries 
such  as  diplomacy,  alliances  and  the  threat  of  direct  intervention.  In  today's  political 
climate,  the  superpowers  are  much  less  likely  to  intervene  in  local  hostilities  with  their 
own  armed  forces;  they  are  much  more  likely  to  attempt  to  support  their  interests  with 
tiie  supply  of  weapons.  Final!},  the  past  decade  has  seen  a  dramatic  increase  in  the 
demand  for  weapons  by  Third  World  nations.  As  new  Third  World  nations  have  been 
created,  and  others  have  continued  to  develop,  they  have  sought  to  acquire  the 
weapons  of  the  industrialized  world.  Given  the  important  role  of  arms  transfers  in 
global  pclifcs.  then,  some  measure  of  control  over  this  activity  would  be  extremely 
desireabie.  .Andrew  Pierre  observes.  ".Arms  transfers  ...  should  be  managed  so  as  to 
prevent  or  contain  confict  and  enhance  the  forces  of  moderation  and  stability." 
.Ref  2:  p.  "I 

One  avenue  of  arms  control  might  be  the  voluntary  restraint  of  the  recipient 
nations.  To  see  if  this  approach  is  plausible,  the  motivation  of  the  recipient  for 
procuring  weapons  must  be  examined.  There  are  a  variety  of  reasons  why  Third  \V‘orld 
v-cuntries  purchase  weapon  systems;  the  need  for  defense  from  external  aggression,  an 
indication  of  regime  legitimacy  and  a  desire  for  world  prestige.  Howeser  it  can  be 
easi!}'  argued  that  regime  survival  is  the  primary  reason  for  buying  arms.  Thus,  it  is 
the  immediate  needs  of  the  Third  World  regimes  that  dominates  their  perceived  needs 
for  arms  imports.  The  regimes  are  generally  not  concerned  with  the  global  problem  of 
e^cui-iting  arm,,  lesels  or  of  the  balance  of  power;  they  simply  want  to  rem.ain  in  power 
and  maintain  their  capabhitv  to  suppress  internal  insurgencies.  In  an  environment 
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such  as  this,  it  is  inconceivabie  that  recipient  nations  would  agree  to  any  form  oi' 
vo!untar>'  arms  import  restraints. 

What  then  of  supplier  restraints?  Might  the  maior  arms  suppliers  agree  to  e\ert 
a  collective  measure  of  control  over  arms  transfers’  The  L'nited  States  and  the  Soviet 
L'nicn  account  for  tlie  majority  of  world  arms  transfers  and  it  is  fair  to  say  that  the 
primary  reas  on  the^e  n.icicns  transfer  weapon  systems  is  for  political  and  ideological 
considerations.  However,  second  tier  suppliers  such  as  France  and  Brazil,  are  heuMlv 
moti'.  uted  by  econonuc  factors.  "France  has  often  been  cited  as  the  best  example  of  a 
country  whose  arms  trade  policies  are  deternuned  more  by  commercial  titan  political 
motivations.  Long  before  the  LD'a  energy  crisis,  the  French  land  to  a  lesser  extent  the 
British!  sought  tlte  most  lucrative  available  markets  and  were  largely  uninhibited  by 
political  restraints.  '  [Ref  1:  p.  6')j  Similarly.  Brazil  exports  95° o  of  its  arms  production 
IRef.  3;  p.  S”].  They  n'liike  no  distinction  regarding  prospective  customers  -  in  fact  they 
currently  soil  arms  to  both  Iran  and  Iraq.  When  the  dichotomy  between  L'S  and 
SoMct  ideclogies  is  considered,  it  is  unlikely  that  these  two  superpowers  will  be  able  to 
agree  on  bilateral  arms  transfer  constraints.  Furthermore,  considering  the  economic 
moti'.  ation  cf  second  tier  nations,  it  is  just  as  unlikely  that  these  countries  would  find 
satisfacter;.'  multilateral  constraints.  Thus,  unilateral  constraint  by  the  United  States 
appears  to  be  the  oniv  likely  avenue  of  exerting  any  measure  of  control  over  the 
transfer  cf  weapons  system.s. 

\\'hat  are  some  examples  of  these  controls  that  the  US  may  attempt  to 
implement?  Gompaf  and  \'ershbow  observe  that.  'It  must  be  recognized  that  some 
tt  pes  of  arms,  by  their  technical  nature,  are  generally  stabilizing,  especially  those  that 
(.an  be  mere  oiTectisely  employed  for  defensive  than  for  clTensive  purposes. "  [Ref  1:  p. 
131  It  may  be  in  the  best  interest  of  the  US  to  promote  the  trans.f'er  of  these  weapons. 
But  v.T.at  actions  might  the  US  take  to  accomplish  this  and  how  elTectivc  would  these 
measures  be?  .Another  exam.ple  cf  a  form  of  control  over  transfer  of  arms  proposed  by 
Hessing,  et  al,  is  to  make  the  transfer  process  more  viscous  by  extending  the  time  I'or 
debate  in  Congress  [Ref  1:  p.  1"<)!.  What  would  be  the  result  of  such  a  policy?  '3\'cu!d 
:t  have  any  effect  on  the  likelihood  cf  translers  to  the  recipient?  To  answer  questions 
like  tl'iese.  the  foreign  policy  analy^t  would  need  some  sort  of  model  of  th.e  arms 
transi'er  process  so  th.at  he  could  test  tiie  various  policies  and  see  what  the  likely  eifcct 
'.'.  .mild  "c. 
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Laurance  and  Mullen  note  that,  'most  of  the  scholars  and  analvsts  conducting 
research  on  arms  transfer  issues  are  generalists  and  political  scientists  draun  to  tiie 
issues  by  concern  for  policy  and  the  policymaking  process.'  [Ref  ?:  p.  S4|  The.  feit 
that  there  is  a  need  for  specialists  from  more  diverse  fields  of  study  to  be  brought  into 
the  investigation  of  arms  transfers.  The  purpose  of  this  study  is  to  bring  Operations 
Resear-'li  tecitniques  to  bear  and  develop  a  model  of  the  arms  transfer  process,  ’['he 
ntcuel  u’lil  be  structured  so  that  it  can  be  used  as  an  arms  transfer  policy  tool  bv 
analysts  with  ver\'  linuted  mathematical  background.  The  technique  chosen  to  model 
arms  transfers  is  cross  impact  analysis.  The  next  chapter  will  introduce  the  concent  of 
cross  impact  analysis  and  KSIM.  an  enhancement  to  the  original  technique.  Chapter  ? 
will  tauor  RSIM  to  model  the  arms  transfer  process  and  will  be  followed  bv  a  chapter 
discussint:  the  results  of  the  s.mulation. 


II.  CROSS  IMPACT  ANALYSIS  REVIEW 


I.NTRODLCTION 

!'u:.;rc'  is  a  disc;pi;nc  that  explores  diilerent  techniques  of  forecasting. 

^a<;c  mission  of  futures  research  is  to  broaden  our  time  horizons  and  enable  us 
li;.  to  anticipate  long  term  cltange  per  se.  but  also  to  see  hew,  by  controlling  such 
es,  ue  ..an  increase  the  range  of  cur  alternatives  and  select  alternati'.'es  likely  to 
.e  a  better  societv  in  both  the  near  and  longer  time  periods."  [Ref  4i  p.  .'<'1 
:s  researc'ners  tend  to  make  a  distinction  between  prediction  and  forecasting.  In 


■.lew.  predwticn  is  an  attempt  to  specifically  identify  the  future  as  a  senes  of 
If  one  takes  this  view  with  the  future  regarded  as  inevitable,  and  noth.ng  could 
le  to  alter  tins  series  cfesents.  then  the  study  and  analysis  of  the  future  tvculd  be 
ess.  Forecasting,  iicwever.  seeks  to  identify  the  alternative  futures,  estiiriate  the 
e  likelihood  of  these  aiternatives  and  investigate  what  controls  exist  to  change 
.eiihcods.  In  order  for  a  forecast  to  be  believable,  it  must  be  supported  by  solid 
icai  methods.  Tiie  doveiepment  of  analytical  tools  for  the  purpose  of  obtaining 
le  I'orecasts  is  the  cornerstone  of  futures  research. 

Most  r'orecasting  teciiniques.  such  as  moving  averages,  regression,  and  '. ariciis 
lung  rcutines.  require  substantial  data  to  implemient.  Howes er.  some  ss  stems  are 
ique  titat  no  hard  trend  data  exists.  Furthermore,  because  the  luture  ci'  some 
IS  IS  so  e\ olutionary.  simple  extrapolation  of  trend  and  growth  curves  is  an 
^fastciw  method  for  developing  a  plausible  forecast.  Sirrulariy,  although  classical 
'lilt;.  ..nd  Statistics  m.ay  deal  with  tune,  it  is  always  in  terms  of  the  past  or 
It  rutiter  than  the  future.  For  this  reason  standard  statistical  techniques  has e 
lound  to  be  lacking  in  their  ability  to  analyze  the  future.  Thus,  in  order  i'or 
'  research  to  be  effective,  it  must  employ  techniques  that  take  adsai.tage  ol 
r.ental  data,  since  most  of  the  i:',k,.Hit  nccessarc  for  the  anticipation  c !  J 1 4 1 X-,  u  c 
-  .'  simply  not  asailable  through  any  ether  means.  .As  Helmer  points  cut,  it 


'e  rcccgni/ed  that  lutures 
f'-  'e  .nmidered  a  part.  ;s  n 


'  minor 


aivsis.  hke  operations  analysis,  of  wliich  it  'hou'.J 
evitable  conducted  in  a  domain  of  u  iiat  might  be 
means  that  dependance  on  intuitive  judgemei.t  n 
t  m  fact  a  mandatory  requirement.'  [Ref  p,  H 
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Futures  research  must  also  involve  inputs  from  fields  related  to  the  area  of  mierest. 
although  not  directly  part  of  that  area.  This  is  because  changes  in  one  area  al'-'C^ien 
may  have  a  "spiil  over"  elTect  and  ,  conversely,  changes  in  other  aspects  of  societ}  m,;y 
alfect  the  area  of  study.  Finally,  the  research  must  be  systematic,  as  a  strict 
methodological  approach  is  often  the  only  tvay  to  sort  tlirough  the  mans  complex 
issues  that  faec  the  futures  researcher.  The  concept  of  using  an  expert,  or  a  contmittec 
of  experts,  in  some  sort  of  structured  analysis  soon  emerges  as  one  the  most  uselul 
methods  available  in  the  held  of  futures  research. 

The  purpose  of  this  chapter  is  to  show  that  cross  impact  analysis  is  a  uselul 
method  of  incorporating  expert  opinion  into  an  analytical  model.  The  chapter  '.'.li; 
present  tite  Delphi  technisTue  and  demonstrate  how  the  original  cross  impact  model  was 
developed  to  rectify  perceived  shortconungs  in  Delphi.  .-X  review  of  cross  impact 
enhancements  will  follow.  This  review  is  not  comprehensive;  it  only  presents  what  this 
researcher  feels  are  the  most  widely  used  enhancements.  Nor  does  the  review  of  the 
various  techniques  go  into  enough  depth  to  stand  alone;  readers  with  farther  interest  in 
any  of  these  developments  are  refered  to  the  original  papers. 

B.  DELPHI 

The  Rand  Corporation  was  prominent  among  the  early  pioneers  in  futures 
research.  They  observed  how  well  expert  opinion  fit  into  the  requirements  described 
above,  but  believed  that  a  group  of  experts  was  eminently  more  capable  of  providing 
an  accurate  forecast  of  the  future  than  a  single  expert.  "The  basic  notion  was  that  the 
group  would  interact  to  com.pcnsate  for  the  biases  of  individual  members,  and  that  the 
hnov,  ledge  of  one  member  of  the  group  may  compensate  for  another's  ignorance." 
iRef.  p.  ISj  Flowever,  assembling  a  group  of  experts  in  one  room  was  not  always 
conducive  to  achieving  the  goal  of  the  group;  namely  arriving  at  the  best  possible 
:'.orccast.  The  reason  was  that  group  dynamucs  often  became  the  prime  motivater  of 
tite  panel,  rather  than  a  common  view  of  the  objective.  Some  of  the  specific  problents 
assoviiated  with  group  dynamics  are: 

1.  .A  group  may  exert  significant  pressure  on  its  members  to  conform  with  the 
group  view,  even  though  the  group  view  may  be  wrong. 

2.  .A  strong  vocal  minority  may  overwhelm  the  majority  by  the  sheer  volume  ol' 
titeir  arcum.ents  wrulc  the  merit  o!' their  arguments,  taken  indw  iduallv.  max  be 
minimal. 


U 
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In  grcups  v^'b.ere  no  on'icial  leader  has  been  appointed,  any  one  indiMdu.i;  nvi’. 
ha'.  e  an  undae  intpaci  on  the  outconte  c:' the  ..ciraruttee  deeiMon  based  or.  ti.e 
st.'ength  of  his  personality,  no:  the  men:  ol'his  arguments. 

Groups  are  often  more  interested  in  reav.r.ing  agreement  on  a  iorecast  rather 
:;n.n  .‘mding  a  t'orecas:  that  may  be  more  accurate.  \et  ollend  some  oomnnttee 
members. 


Ire  groan',  as  a  uho! 
eo..^,.:.  m,  and  hel^ 


;  '.vii;  prcbaniy  snare  a  Lcmmcn  nais  based  on  its  eiiiture, 
:f  expert. se 


I'ne  Delphi  prooedare  of  foreeasting  teas  deselcpeJ  by  the  Rand  Corporation  as 
.1  may  fo  eammate  most  ot'  :;ie  disadv antages  ot'  group  dynanuos  mhiie  eapit.di/ir.g  on 
'.'me  v-ohcLtis  0  kno.v'.edge  inherent  in  a  group  of  experts.  T'ne  entire  Delpiu  procedure 
■.'..;s  ^or.du..ted  '.la  a  senes  of  siuesnonnaires;  the  panel  veas  ne\er  reoaired  to  be 
asscn.-'led  m  a  ^.cmmcit  looation  for  tine  procedure  to  u'ork,  d'he  identity  of  group 
me;r.''ers  was  not  made  hno’.vn  to  individual  members  which  avoided  a  spesiik  opinion 
'•'emg  linked  to  a  specific  panelist.  The  questionnaires  were  structured  in  order  to 
pro  . '.do  tecdback  to  the  group  at  each  step  of  the  procedure.  .\t  the  conclusion  of  the 
Deir'm,  t;:e  director  had  a  statistical  sun'umary  which  included  the  opinions  of  the  entire 


.•\  brief  outiin 


e  c; 


rh.s 


issic’  Delphi  as  described  by  Martino  (Ref  t.'!  follows: 


Step  1.  Par.ehsts  are  asked  to  provide  a  forecast  of  the  future  in  the  subject  area  of 
interest.  Tins  s;ep  ;<  completely  unstructured;  some  panelists  may  provide  a  narratise 
Scenario  wlnie  others  may  provide  a  list  of  events  and  associated  dates.  The  director 
must  ta.xe  these  inputs  and  consolidate  them  into  a  final  list  of  key  events  that  are  as 
spec.fic  as  possible. 

.Step  2.  l  i'.is  ;i  st  of  events  is  returned  to  the  group  who  are  asked  to  estimate  the  date 
by  which  the  event  will  have  occured.  In  addition,  the  group  members  must  state  a 
reason  whv  they  chose  that  date  as  their  estim.ate.  The  director  then  prepares  a 
stmist.cal  sunmiary  of  the  groups  opinions  (miedian  and  inter-quartile  range  for  each 
oven:  I  and  summarizes  the  arguments  supporting  the  date  estimates. 

Step  The  Delphi  panel  is  new  presented  with  a  list  of  events,  the  statistical 
■'Umm.iry.  and  reasons  for  date  estimates.  Individuals  are  asked  to  renew  ihe 
.u'gumei'.ts  and  make  a  new  esi.m.ne  for  the  date  r\  '.vhich  an  event  w.il  liase  o^i.'..rred- 
1:  ’he  r.^.v  es:;m.;te  tails  outside  the  mter-quarthe  range,  the  panelist  :s  asked  to  .'U'-tifv 
■■:cv.  .i:;d  ..ur.men:  on  opp 's;ng  arguments  Ihe  u;rcv.t.m  no.',  ^ampmcs  ne.'. 
...eh.  .r.'  .r.h  a  w.le  riners  ■  :iie  urh.i:ch  csurnates  .aid  .icun  ^..n.n.ari.'cs  y-.e 


timJimiliidCmJi 


Step  4.  The  updated  medians  and  inter-quartile  ranges  are  presented  to  the  panel  for 
one  more  evaluation.  .\s  in  Step  3.  if  estimates  fall  outside  the  inter-quartile  range, 
panelists  must  justify  their  extreme  vic.vs.  The  director  can  now  compute  a  final 
median  and  inter-quartile  range,  and  he  has  a  set  of  arguments  germane  to  these  events 
;cr  which  no  date  was  settled  on.  The  forecast,  then,  consists  of  the  list  of  events  and 
tneir  median  and  inter-quartile  dates. 

O.fen.  there  were  events  included  in  a  Delphi  analysis  that  were  inter-reiated. 
T’ne  v.ay  tiiat  the  procedure  was  structured  required  that  estimates  for  the  occurrence 
of  one  of  these  inter-relatcd  esents  be  given  without  consideration  for  the  dates  of  any 
other  esents.  If  the  occurrence  of  one  event  was  predicated  upon  the  prior  occurrence 
cf  another  event,  this  proved  to  be  a  serious  problem  in  providing  accurate  estimates. 
It  was  I'or  this  specific  reason  that  cross  impact  analysis  procedure  was  created.  "The 
genesis  of  cross  impact  was  the  problem  that  Delphi  panelists  were  sometimes  asked  to 
make  i'orecasts  for  indi\idual  events,  when  other  esents  in  the  same  Delphi  could  affect 
tiiese  e'.ents.  Thus,  it  was  recognized  that  there  was  a  need  to  allow  for  these  cross 
impacts  :rom  one  event  to  another."  [Ref  p.  61] 

C.  CROSS-IMP.ACT  A.NALVSIS 

Cross  impact  analysis  was  first  developed  by  T.  J.  Gordon  and  H.  Hav.vard  in 
f'lhS,  Their  goal  was  to  overcome  the  difficulty  of  event  interrelationships  by 
estimating  and  explicitly  accounting  for  them.  They  reasoned,  "Most  developments  are 
in  some  way  connected  with  other  events  and  developments.  It  is  hard  to  imagine  an 
event  witlicut  a  predecessor  that  made  it  more  or  less  likely  or  influenced  its  form  -  or 
one  which,  after  occurring,  left  no  mark.  This  interrelationship  between  events  is  called 
^rcss  i.mpact".'"  [Ref  8:  pp.  100-101]  Events  connected  by  cross  impacts  had  two 
modes;  even:  .A  couid  either  enhance  or  inhibit  the  occurrence  of  event  B.  These  event 
linkages  were  then  assigned  a  strength  such  that  a  high  strength  indicated  the 
occurrence  of  event  .A  exerted  a  large  inllucnce  on  the  probability  of  occurrence  of 
even:  B.  Using  this  data,  a  cross  impact  matrix.  Sy.  was  developed  showing  how  the 
occurrence  of  every  event  j  impacted  on  the  probability  of  occurrence  of  every  other 
ev  en:  i.  Once  an  initial  probability  of  occurrence.  P(i),  for  each  event  was  estimated, 
the  mode!  was  ready  to  be  tested. 

The  cro'-s  immact  mode!  developed  by  Gordon  and  Havvvard  was  stochastic  in 
nature,  fhc'.  felt  that  after  a  laree  number  of  runs  (on  the  order  of  1000).  steadv  state 


probabilities  for  each  of  the  events  would  be  reached.  The  basic  procedure  is  described 
here. 

Step  1.  Chose  an  event,  E^.  at  random  from  the  list  of  events.  Using  a  uniform 
random  number  generator  (0,1 1,  determine  whether  or  not  E-  occurred  based  on  its 
mitia!  probability  of  occurrence  P(i). 

Step  2.  If  E,  did  not  occur,  discard  it  from  the  list  and  choose  another  event. 
Continue  this  process  until  an  event  is  found  to  have  occured. 

Step  3.  Modify  the  probabilities  of  occurrence  of  all  remaining  events  according  to  the 
equation: 

P'i  1)  =  P(j)  -  P(j)(l  -  P(j))Sij 

Step  4.  Return  to  Step  1  as  many  times  as  necessary  until  no  events  remain. 

This  four  step  procedure  constituted  a  single  run  and  was  repeated  1000  times.  The 
final  probability  of  occurrence  for  each  event  was  estimated  by  noting  how  many  times 
an  event  occured  in  the  1000  runs  and  dividing  by  1000. 

Gordon  and  Ha>-i.vard  then  tabulated  the  results  by  event,  initial  and  final 
probabilities,  and  a  ranking  according  to  initial  probability  and  change  in  probability. 
From  these  tabulations,  they  drew  conclusions  about  which  events  were  most 
susceptible  to  change  and  which  were  most  resistant  to  change.  "The  ranking  by 
probability  shift  is.  in  essence,  a  list  of  the  items  most  alTected  by  the  suspected 
interaction.  In  other  words,  the  item  which  had  the  highest  probability  shift  could  be 
expected  to  be  the  one  most  influenced  by  external  events  depicted  by  the  remainder  of 
the  list."  [Ref  S;  p.  108]  One  of  the  real  values  of  cross  impact,  though,  was  found  to 
be  the  ability  to  test  the  effect  of  various  policy  decisions  on  the  final  occurrence 
probability  of  the  events.  The  evaluation  of  policy  decisions  could  be  implemented  by 
vaiying  the  probability  of  one  or  more  events,  replaying  the  matrix,  and  comparing 
these  results  to  the  original  results. 

D.  CROSS  IMPACT  IMPROVEMENTS 
1.  Enzer 

While  cross  impact  had  a  great  deal  of  intuitive  appeal,  there  were  significant 
drawbacks  to  the  original  approach.  Gordon  and  Hayward  noted.  "We  beliese  that 
this  work  is  only  indicative  of  a  methodology  of  cross  impacts.  If  possible,  its  current 


shortcorfungs  should  be  corrected  in  future  work.."  [Ref  8:  p.  115]  Sehcwn  Enzer 
recognized  some  of  these  shortcomings  and  sought  to  correct  them.  Enzer  observed 
that  since  the  relationship  between  an  event's  initial  probability  and  its  final  probability 
was  based  on  a  quadratic  f'unction.  it  was  impossible  for  an  inhibiting  event  to 
significantly  change  a  high  probability  and  an  enhancing  event  to  significantly  change  a 
lo'v  probability  event.  Moreover,  the  etTect  of  opposite  cross  impacts  were  not 
symmetric.  For  these  reasons,  Enzer  abandoned  the  quadratic  manipulation  suggested 
by  Gordon  and  Ma^-ward  and  developed  a  cross  impact  method  based  on  the  likelihood 
ratio.  IRe.f  4]  .An  o’.  erview  of  this  model  follows. 

The  odds  of  an  event  occuring,  Oii),  are  compvued  from  the  probability  of  occurrence. 
Pi  n.  as  Ibllows: 

0(1)  =  P(i)  [1  -  P(i)] 

Similarly,  probability  can  be  computed  from  the  odds  by 

P(il  =  0(1)  [1  ^  0(i)] 

If  event  j  occurs  and  changes  the  odds  of  event  i,  then  the  likelihood  ratio,  R-  is 
related  to  the  odds  of  i  given  j  occured  by  the  following: 

0(i,j)  =  R-O(i) 

This  implies  that 

P(i!j)  =  0(i|j)  (1  +  O(ili)] 

Using  these  relationships,  then,  the  modified  probability  of  the  occurrence  of  i,  given  j 
has  occured,  is 

P(i])  =  RyP(i)[l  -r  (Rjj  -  l)P(ii]  (eqn2.li 

Enzer  used  the  same  basic  .Monte  Carlo  approach  of  Gordon  and  Ha\-ward. 
but  modiiied  the  probabilities  after  cross  impact  by  equation  2.1  Through  the  use  of 
the  likelihood  multiplier,  the  magnitude  of  change  in  initial  probability  is  the  same  for 
reciprccai  likelihood  ratios  and  'the  domain  of  change  that  is  permitted  by  this 
tec'nnique  extends  from  0  to  infinity,  so  that  impacting  events  can  have  the  elfect  of 
totally  eliminating  the  possiblity  of  occurrence  of  a  subsequent  event,  or  in  lact  causing 
tite  subsequent  event  to  occur.  "  [Ref  4:  pp.  43-44]  Through  the  use  of  a  likelihood 
r.itio.  En/cr  '.vas  abie  to  overcome  th.e  problems  of  assymetry  and  the  diminishing  elTect 
m  tltc  earhewt  form  olGrcss  ur.nact  anah. s;s. 


Turoff 


Murray  TurotT  approached  the  same  set  of  discrepancies  in  much  the  ^ame 
way  as  Enzer.  except  that  TurofT  defined  an  occurrence  ratio.  0),.  as  the  natural 
logarithm  of  the  odds.  [Ref.  9j  That  is. 

0),  =  (DfP.)  =  In  O;  =  InIP,  1 1  -  P,)]  (eqn  2.:i 

Lsing  this  deiinition.  TurolT  set  about  establishing  a  relationship  between  the 
likeliiiood  of  occurrence  of  an  event  and  the  elTort  put  forth  to  enhance  or  inhibit  that 
coourrence.  He  wanted  this  relationship  to  be  "such  that  if  an  equal  amount  of  eiTort 
is  devoted  to  both  enhancing  and  preventing  the  occurrence  of  an  event  then  the 
likelihood  corresponds  to  a  probability  of  one-half  (i.e..  random  or  neutral)." 
:Ref.  "V-  p.  319]  Furthermore.  Turoff  assum.ed  that  the  estimater  would  provide 
^or.sistent  probability  estimates;  this  was  a  crucial  dif'erence  between  his  approach  and 
tl'iOse  techniques  using  Monte  Carlo  simulations.  With  these  assu.mptions.  he 
com.puted  an  expression  for  P^.  the  probability  of  occurrence  of  event  i: 

I 

P;  =  - = -  leqnd..' 

1  m  e.xpl-v.  -  lqj,Pi,) 

k^i 

where  y-  was  a  function  of  unknown  variables  and  C|j.  was  the  cross  impact  term. 
TurolT  then  showed  that  the  occurrence  ratio  could  be  related  to  equation  2.3  by  using 
equation  2.2  to  get 

*  Yi  -  IC.kPt 

He  noted  that  while  the  impact  of  the  k^^  event  on  the  i^*^  event  was  additive  with 
respect  to  the  occurrence  ratio,  it  was  multiplicative  with  respect  to  odds: 

Oi-Pj  ‘  AP^)  =  0;aPj)0,j.APp 

In  o'i'.er  ’.vords,  any  change  in  the  probability  of  occurrence  of  event  j  which  had  an 
cdcw;  on  c. ent  ;  would  change  the  odds  multiplica:i\ el>'.  He  concluded  his 
dc. eiopnicn:  '.sith  the  notion  that  his  series  of  equations  satisfy  a  likelihood  ratio 
•  .c’.'.  point  o;'  statistical  inference  since  the  final  odds  of  an  esent  occuring  may  be 
vr.tten  as  ti.e  product  of  an  initial  odds  times  a  likelihood  ratio. 
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3.  Helmer 

Enzer  and  TurofT  utilize  transformations  into  odds  space  before  collectiseis 
applying  the  cross  impacts.  The  inverse  transformations  alwass  map  onto  the  closed 
interval  from  0  to  1.  thus  making  it  impossible  to  create  an  illegal  probability.  Olaf 
Helmer  proposed  transforming  probabilities.  P.  into  'R-space  "  to  accomplish  th.e  same 
Sanction  c:  mapping  <0.1)  onto  itsell  in  his  modification  to  the  cross  impact  method. 
;Ref  0] 

Helmer  noted  that  trends  have  a  natural  lower  bound  at  0  and  that  the 
ma.Mmurn  natural  upper  bound  of  all  the  trends  could  be  used  as  the  upper  bound.  L  . 
lor  a.i  trenus.  He  a. so  needed  a  central  value”,  C,  which  he  let  equal  the  median  of 
t.ie  trend  s  estimated  value.  With  these  parameters.  Helmer  defined  his  transformation 
as 

P  -  C 

R(  P )  =  R -  ,  0  <  P  <  L” 

PfU  -  P) 

where  K  -.vas  defined  as 


K  = 


(C  -  S)(L'  -  C  -  S) 


s 


J.  Pi 

■.  p.t 


d  S  '.vas  a  'surprise  threshold".  The  surprise  threshold  was  set  such  that  "we  will  not 
surprised  depending  on  whether  the  true  value  turns  out  to  be  outside  or  inside  the 
er-cuartile  interval. '  (Ref  5:  p.  23]  In  this  model,  cross  impacts  are  additive  in  R- 


space: 


RfPfiij))  =  R(P)  +  Xij  on  i) 


svnere  .Xij  on  ii  ;s  the  actual  cross  impact  coefficient  from  the  cross  impact  matrix.. 
Helmer  recommended  a  scale  of -3  to  ^3  for  the  values  of  X  according  to: 

±12  =  small  ±  2  =  large 

±  1  =  medium  ±  3  =  overwhelming 

Po'itive  and  negative  values  corresponded  to  enhancing  or  inhibiting  cross  imrsacts  as 
before.  Once  all  impacts  had  been  aggregated,  one  need  only  apply  an  inverse 
transformation 


I  RE  -  K)  *  v'RL  -  Ki-  ^  4RC 


to  obtair.  the  finai  probabiiitie';-  A  Moriihcar.t  ditTcrence  in  Fielrner  s  r.'.cdei  and  titC'-c 
preseiUeJ  earlier  :s  that  ;Ir.pa..t^  la't  lor  onb.  one  time  interval,  but  this  \>.as  oorre^'cd 
the  Use  oh  a  carr.  o.cr  paran.eter  lOiivn  iielmer  ue'-i.ribcd  in  [Ref.  'J. 
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if 
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er  doe',  a  <.ipe;-b  joA 
10  O',  aiuat.on  of  'i.c 
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'pear-,  .n  jReh  l>t|.  ir.  '.hiv  work.  A1 
o:  'iiiii;. 0  'iie  .por.  a..:iC'  aiiu  oonduLtin'O  a  detaiiod  anal'. 

1.0.:  Motiio.;'  lie  e.aiu.ite'.  'he  n.odeS  ''a'od  on  tiieir  interna;  ooi 
renlien..;  .  :r,;e;h:;.,  oener.dit;. .  .ind  daritv.  Alters  paper  should  deliniteiy  be  consulted 
i;  there  >  :..rtlter  inierc'.t  in  eiti.er  the  likelihood  rat.o  model  or  the  R-^pace  mode: 
d  Bloom 


.A  more  ^traiehtlorward  but  somewhat  restrictive  cross  impact  model  was 
d  Ay  MitvOitei  Bloom  in  R)~5.  [Retl  11)  In  this  model,  the  objectite  was  to  u^e 
nt  past  to  loremast  the  near-term  future  of  the  system  under  consideration. 
1;. ,  a  prerequisite  for  the  use  of  Bloom  s  model  is  the  availability  of 
a'l'ie  data  and  a  system  for  study  which  lends  itself  to  data  collection.  There 
'■'as  asstintptions  upon  which  the  model  is  built: 


The  system  possesses  considerable  inertia.  The  most  likely  future  of  the  system 
'.vhl  be  a  simple  extrapolation  of  past  variable  trends,  ^^’hlle  the  probability  of 
this  ocvurrm.g  may  be  small,  it  is  assumed  to  be  larger  than  anv  other 
.utcrnative  luturc. 

fne  deviation  of  any  trend  from  its  undisturbed  future  path  will  hase  a  cross 
impact  on  th.e  other  trend  \ariables  m  the  system. 


fitere  are  several  reasons  why  Bloom  felt  that  a  deterministic  trend  cross 
im.pa^.:  mode!  with  a  graphical  display  of  system  \ariables  was  more  appropriate  than 
me  models  presented  earlier,  f  irst,  cross  impact  models  produce  a  I'lnal  set  of  esent 
proAubiiitiC',  they  do  not  yield  a  likely  scenario  that  resulted  in  those  probabilities,  nor 
do  ti'.cy  provide  the  user  with  the  dynamic  change  in  the  trend  variables  oser  time. 
Second.  Bloom  noted  that,  "the  cross  impact  method  has  not  dealt  with  system  trends, 
..ilthcugh  It  has  been  asserted  by  students  o!'  social  change  that  persistence  'Aitiun 
'Cs.et.ii  systems,  as  indicated  by  trends,  is  equally  important  in  understanding  h>tor. 
-.s  st.dden  change  caused  by  singular  e.ents,”  (Ref  11:  p.  3S]  F’lna'.iy,  a  trend  sross 


.mp.n.:  m.;dei  is  a  way  to  combine  past  data  with  the  intuitive  perceptions  cl  a  group 
.  :  expert',  on^e  .i  graphical  picture  of  the  future  is  often  easier  for  a  group  to  emihiate 
',!  :>  I'.olistie  mow  of  the  future.  \S'ith  this  rationale.  Bloom  dc'- eloped  an 


Vi 


equaf.or.  that  is  easily  used  by  a  panel  of  experts  with  hrruted  background  in  '.he 
mathematics  of  forecasting  and  who  have  access  to  relevant  data. 

Bloom  s  basic  cross  impact  ecu.iticn  is 

X'(t  -  At!  =  X,it)  -  AXg  1  * 

1 

where  Xht)  is  the  level  of  variable  X:  at  time  t,  AX,  is  the  extrapolated  chanee  in  X, 

1  ]  "  j 

over  it.  t  -  At)  and  g-  is  the  cross  impact  of  trend  j  on  trend  i.  Equation  2.4  is  mere 
easily  understood  if  written  in  the  formi 

Cross  impacted  value  of  a  prion  change  due  to  the 

value  of  trend  =  trend]  increment  x  cross  impacts  of 

1  at  It  -  At  I  at  time  t  to  trend  ]  all  other  trends 

There  are  several  points  to  be  made  about  Bloom's  approach.  Equation  2.4 
expresses  the  future  of  a  system  in  terms  of  persistence,  not  in  terms  of  change.  "It 
explicitly  excludes  crises  and  other  system  shocks  which  result  in  sudden,  large, 
discontinuous  changes  in  the  levels  of  the  key  variables  in  the  system.  "  (Ref  II:  p.  .''2| 
^'aria''lcs  must  be  vie'wed  as  those  which  change  monotonicly.  The  model  cannot  deal 
'.Mth  trends  tiiat  are  constant  ever  time  since  this  will  make  AXj  equal  to  zero  and 
implies  that  cross  impacts  have  on  effect.  Finally,  [Ref  11)  is  restricted  to  only  trend- 
t.'cnd  impacts.  This  last  deficiency  was  corrected  when  Bloom  developed  an  extension 
to  hi'  original  model  to  include  events  [Ref  12).  The  inclusion  of  events  is 
a..c;ontp!:'hed  by  using  the  cumulative  distribution  function  iCDF)  for  the  logistics 
c.'tri'i'aticn  to  describe  the  probability  trend  of  an  event  over  time.  Consequently. 

ent'  cct.ld  be  viewed  as  trends  tinat  increased  monotonicly  from  0  to  som,e  maximum 
..hue  O’. er  a  range  of  time  determined  bv  the  parameters  of  the  equation.  While  it  is 
rec.ocm/ed  that  a  true  CDF  must  approach  a  value  of  one  as  t  approaches  infinity,  the 
-M-a'ion  needed  is  obtained  bv  miultiplvmg  the  logistics  CDF  by  a  maximum  value  \1. 
r  SMIC-74 

.\  vtoss  impact  model  that  is  substantiallv  dilTerent  in  approach  and  results  is 
SMK'-'-  ^h  was  ^cncer-ed  bv  Dupernn  and  Godet  in  19''5.  [Ref  1.']  Thev 
■■.'r-eu  :m.>;  o!  tl'.e  .rnpam  methods  developed  up  to  this  point  failed  to 


rc'. '.dc  consistent  final  nrobahihties  and  vieMed  results  such  as 


Pt  1)  Pi  i  ]  iP'  j  ■ 


■;'.ic:;  is  incomnatihle  with 


P(i)  =  Pii  j'Pl;)  -  Pm  net  jiPinet  ]\. 

)upcrrin  and  Godet  noted  that  for  a  system  of  n  events,  there  were  2’^  centhinaticns 
.  c. ci'.ts  possible.  1  his  could  be  viewed  as  2'^  scenarios  whw'n  ccrrc’spcndea  to  e\e:',:  ; 
cc..rr:i'.g  cr  not  occurring  for  all  events  from  one  to  n.  Th.e  most  prerabie  '.Ci'.ario 
sua'.iy  i;ad  a  prchability  of  occurrence  of  only  about  "  I,  depending  on  the  nuiv..Mr  cf 
■vents  .;nd  the  initial  probabilities  assigned  to  the  events  iP>.ef  1':  p.  eOdf  However, 
'.ere  were  usually  vpuite  dillerent  scenarios  with  only  slightly  lower  probabilities  tiuin 
te  incst  probable  scenario.  Duperrin  and  Godet  felt  that  for  cross  impact  to 
verconte  this  drawback,  the  method  employed  would  have  to  produce  a  rank  order  cf 
il  possible  scenarios.  fhe  requirement  to  generate  a  consistent  set  of  final 
robabiiities  and  to  produce  a  scenario  ranking  led  to  the  development  of  SMlG-'4. 

The  method  assumes  that  a  panel  of  experts  will  be  able  to  render  opinions 

^out; 

•  fite  list  of  n  events  which  are  considered  to  be  the  key  ones  for  the  system 
under  consideration. 

•  The  probability  P(i)  of  each  event  defined  as  the  probability  of  the  occurrence 
of  Cj  within  the  time  period  considered. 

•  The  conditional  probability  of  the  separate  event  pairs: 

Pm']!  =  the  probability  of  i  if]  occurs 

Pm  not  ])  =  the  probability  of  i  ifj  does  not  occur 

.A  given  scenario  denoted  Ej^  is  composed  of  n  separate  events;  or  ebar;  for 
d  events,  which  correspond  to  event  i  occuring  or  not  occuring,  respectively.  This 
^enario  has  an  unknown  probability  Kk  and  the  sum  of  ail  must  equal  one.  P'  m'  is 
.eiir.ed  as 

PH  11  = 


^ik  ~  *  ’■*  “^i  ^k 

6;^  =  I)  if  ebar  forms  part  of  Er 


['or  more  uetails  on  tlie  issue  cf  consistency  in  cross  intact  analysis,  sec  [Ref  14] 
nd?  Ref  151 


.'-•V-'.'.'.  a- 
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Similar  definitions  are  established  for  P'd'i  j)  and  P*(i']ban.  These  theoretical  values  of 
P'vi).  P*(i:j)  and  P*iTjbar)  must  satisfy  the  constraints  established  by  Baysean 
probability  theoip-.  Duperrin  and  Godet  wrote  an  objective  function  which  was  the 
diiference  between  the  Piij)  estimates  from  the  experts  and  the  theoretical  i  j  iP  ■'(] ) 
factors  expressed  in  terms  of  rr|^.  This  objective  was  mimnuzed  subject  to  the 
constraints  titat: 

^Tt^  =  I  and  JT^  >  0  for  all  k 

The  minimization  problem  described  above  is  of  the  quadratic  form  with  linear 
constraints.  Solving  this  svstcm  of  equations  yields  both  a  consistent  set  of 
probabilities  and  a  cardinal  ranking  of  all  possible  scenarios. 

S.MlC-~4  came  under  quick  criticism*  by  Mitcliell  and  Tydeman  [Ref  16].  They 
were  able  to  show  that  the  cardinal  ranking  of  scenarios  was  not  unique  and  further, 
that  the  system  of  equations  is  potentially  very  large.  If  the  problem  was  reformulated 
in  terms  of  P’Ni»  and  P*!]).  and  appropriate  constraints  added,  then  the  system  would 
becom.e  a  linear  programming  problem  which  could  identify  the  multiple  solutions. 
Furthermore,  't'.  Ka>a.  et  al  (Ref  I":  p.  245|  show  that  the  use  of  Dupperin  and 
Godot's  quadratic  objective  requires  such  computational  effort  that  obtaining  a 
solution  IS  prohibitive  if  n  is  very  large.  Finally,  the  ability  of  expert  panelists  to 
answer  questions  that  lead  to  estimates  of  P(i:j)  and  P(ii]bar)  was  doubtful.  .Mitchell 
and  Tydeman  pointed  out  that  "results  of  studies  currently  in  progress  indicate  that  ( 1 1 
participants  arc  frequently  confused  and  unsure  of  the  interpretation  of  such  questions 
-ind  '2)  respondents  often  interpret  the  questions  in  terms  of  time-dependant 
*;ondi:ionai  probability  statements."  (Ref  18i  p.  133]  The  significant  difference  between 
P' 1  j  I  and  Pdj  occured  first)  is  clearly  demonstrated  in  (Ref  18).  Consequently, 
SMIC-"4.  which  appears  to  be  a  fairly  popular  cross  impact  approach  has  some  major 
diiTiculties  '.vhich  hamper  its  effective  use. 
fx  KSIM 

In  1972,  Julius  Kane  advanced  another  cross  impact  technique  that  was 
similar  to  several  of  the  techniques  already  presented,  yet  was  unique  in  its  approach 
,Ref  19j.  Unlike  other  futures  researchers,  with  the  exception  of  Bloom.  Kane's 
primarv  concern  was  to  develop  a  model  that  could  be  used  by  people  with  little 
matb.ematicai  background.  Me  believed  that  most  simulation  models  were  e.xcessiveh 
numerical  and  focused  attention  on  those  variables  which  were  easily  quantifiable  and 
tended  to  exclude  '.ariables  that  were  basicly  subjective  in  nature.  .-\s  a  result,  mo'^t 
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policy  niakcr';.  lor  uiuuu  ihc<;c  inodch  wore  dcNigircd,  were  rchittani  lo  invest  the  time 
to  uiidci'statid  how  to  rise  the  <;iimilatioii  model  Kane  s  objective  was  to  attempt  to 
rcctily  this  situation.  "It  was  the  purpose  of  our  research  to  try  and  design  a 
simulation  procedure  -  or  better  yet.  a  simulation  language  m  which  teclinically 
unsophisticated  people  could  quickly  become  fluent  in  the  logical  expressions  of  cross 
impact  concepts.  "  [Ref  I'l:  p.  lyo)  Additionally,  Kane  sought  to  structure  the  problem 
to  allow  for  a  realistic,  graphic  display  of  the  sc  stem  variables  since  this  was  most 
easily  processed  by  the  user.  'Ihese  notions  form  tlie  basis  of  KSI.M,  Kanes 
SIMulation  language,  fhe  model  is  constructed  such  that  it  has  the  following 
properties: 

•  System  variables  are  bounded.  With  an  appropriate  set  of  units  these  bounds 
can  always  be  set  at  u  and  l.U. 

•  .A  variable  increases  or  decreases  according  to  whether  the  net  impact  by  other 
system  variables  is  positive  or  negative. 

•  I  lie  response  of  a  variable  to  a  given  impact  decreases  to  0  as  that  variable 
approaches  its  upper  or  lower  bound. 

•  A  variable  will  produce  a  greater  impact  on  the  system  as  it  grows  larger 

•  Complex  interactions  can  be  broken  down  into  a  network  of  discrete  binary 
interactions. 

The  KSIM  algorithm  uses  state  variables.  X^ft),  which  are  bounded  by  0  and  1.0  as 
desenbed  in  the  lust  model  property.  The  updated  value  of  Xj(i)  is  calculated  by 

^i(t) 

Xj(t  +  At)  =  Xj(t)  (eqn  2.5) 

where  Tt^ft)  is  chosen  explicitly  as 

I  +  '■2AtX(!Y,j(t,)|  -  Yjjit)) 

j 

K|(t)  =  -  (eqn  2.6) 

1  +  CAtXdYijitll  ^  Yijlt)) 

J 

and  Y,j*T)  arc  the  cross  impacts  of  variable  j  on  variable  i  and  At  is  the  tunc  period  of 
one  iteration.  The  cross  impacts  arc  functions  of  the  magnitude  of  tlie  impacting 
\jr::ibie  and  the  rate  of  change  of  that  variable  as  shown  in  equation  2.7. 


•t  f 


y.  (t)  =  A..X.(t>  ^ 

'll  ■)  j 


B..(dx.(t)  dn 

1)  j 


I  eqn  2.") 


Equation  2.6  implies  that  JTat)  0.  thus  equation  2.5  always  maps  the  open  interval 
1,0.1}  onto  itself.  Equation  2.6  is  much  clearer  if  written  as 


1  -  AtISL'M  OF  NEG.ATIVE  IMP.-XCTS  ON  X;| 

'  '  1  -  AtlSLM  OF  POSITIVE  IMP.ACTS  ON  X.) 

When  the  positive  impacts  are  greater  that  the  negative  impacts,  then  ;t;(t.)  <  1  and 
results  in  an  increase  of  Xa.t).  Similarly,  when  the  negative  impacts  outweigh  the 
positive  ones,  Xjit)  will  decrease.  This  is  in  accordance  with  the  second  assumption 
property  of  the  KSI.M  model.  To  investigate  how  the  remaining  properties  are 
satisfied,  let  At  -»  0.  Equations  2.5  and  2.6  can  be  written  as  a  limiting  system  of 
di.Terential  equations: 


dXj.t\ 

dl 


j 


Kane  refers  to  the  term  Xj(t)lnXj(t)  as  the  "modulator"  since  as  Xj(t)  -*  1.  dXl^i)  di  -* 
h  and  likewise,  as  X|(t)  -♦  1.  InXjit)  -♦  0  and  dXJ^n  di  -*  0.  This  is  the  requirement  of 
the  third  property.  Finally,  if  the  effect  of  Xj(t)  in  equation  2.7  is  considered 
individually,  it  can  be  seen  that  it  will  have  a  greater  effect  on  7tj(t),  and  ultimately  on 
X.'t  -  At),  as  the  magnitude  ofXj(t)  increases.  Thus,  the  fourth  property  is  satisfed. 
Finally,  since  the  system  is  modelled  through  the  cross  impact  term.  y...  which 
describes  a  binary  interaction  between  state  variables  i  and  j.  it  can  be  seen  that  the 
final  KSIM  property  holds.  The  output  of  KSIM  is  a  plot  of  each  of  the  X-  variables 
over  time  and  provides  the  most  concise,  yet  descriptive  means  of  conveying  the  results 
ol’  the  simul.uion. 


Lipinski  and  Tydeman  proposed  an  extension  to  Kane's  model  which  allows 
for  the  inclusion  of  events  as  state  variables  [Ref  20).  Their  extension  is  similar  to 
Bloom  s.  in  that  events  are  described  by  thcir  CDF.  Since  trends  in  KSIM  are 
described  by  logistics  curves,  it  seems  natural  that  a  logistic  CDF  would  ft  nicely  into 
the  Iramework  of  the  model.  Lipinski  and  Tydeman  show  that  if  X^  is  a  trend  variable 
which  •.‘i  initially  ncnconstant.  then  updated  values  of  this  variable  are  calulated  as 


X,(t  ‘ 


At) 


=  X,(t 


,Qj(t)7t,(tl 


24 


'■•C-'v  . 


where  TT/ti  is  the  cross  ir.'.pact  term  ucscribed  previously  and  Q-(t)  is  a  function  that 
generates  the  initial  nonccnstant  trendi 

X  It  -  At)  =  X  itt^'-^--^^' 

fhus,  to  include  an  event  in  the  simulation,  the  function  Q(t)  which  describes  that 
c".  e.".:  St  [i)F  is  Jaund  and  included  in  the  calculation  at  the  appropriate  time. 

E.  E\ALL.\TION 

I  has  researcher  believes  that  KSIM.  in  its  extended  form,  is  the  m*ost  com.plete 
and  inttiitr.eiv  pleasing  model  available  from  the  cross  impact  techniques  reviewed.  It 
siescribes  th.e  complex  interaction  of  variables  in  a  realistic,  nonlinear  fashion,  and. 
once  programmed,  does  not  require  mathematical  sophistication  to  manipulate  the 
various  parameters  and  interpret  the  results.  The  growth  characteristics  assumed  can 
be  mcd.fied  to  accommodate  event  CDF's  of  any  distribution  so  the  model  is  extremeiv 
f.e\f''le.  The  nature  of  the  model  encourages  investigating  the  implications  of  dilTerent 
pc:;c>  dcv-isions.  .Additionally  the  group  interaction  that  results  from  evaluating  inputs 
a.nd  dis,.u<'ing  outputs  is  often  just  as  valuable  as  the  actual  results  of  the  simulation, 
f  or  these  reasons.  KSIM  appears  to  be  a  useful  tool  for  trying  to  find  answers  to  some 
of  ti'.e  questions  raised  in  chapter  one,  and  investigate  issues  surrounding  arms 
trans.'ers.  The  next  chapter  will  develop,  in  detail,  a  KSIM  miodel  that  can  be  used  to 
s'.m.uiate  arms  transfers. 
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III.  DEVELOPMENT  OF  AN  ARMS  TRANSFER  MODEL 


A.  INTRODLCTION 

The  prcMJUs  er.apter  introduced  the  field  of  futures  research  and  described  the 
forecasting  technique  of  cross  impact  anahsis.  Several  imorovements  to  the  basic 
model  '.vere  discussed,  among  them  KSIM.  KSI.VI  is  unique  in  that  it  is  a  new 
mathematical  language  that  facilitates  a  non-technical  decision  maker's  active 
invoivemierit  in  modelling  system  behavior.  This  is  critical  for  investigating  policy 
implications  in  the  area  of  arms  transfers  since  there  is  very^  little  hard  data  with  which 
one  may  construct  a  simulation  model.  The  purpose  of  this  chapter  is  to  develop  an 
interactive  policv  sirniulation  based  on  KSIM  to  describe  the  arms  transfer  process. 

B.  K.SIM  PHILOSOPHY 

Before  proceeding,  it  is  necessary  to  discuss  the  underlying  philosophy  of  the 
snmulation.  The  output  of  a  KSIM  simulation  is  a  graphical  display  of  system 
variables  showing  how  they  change  over  time.  As  such.  KSl.Vl  emphasizes  the 
geometry  of  relationships  rather  than  hard  numerical  predictions.  As  an  e.xample  of 
what  is  meant  by  the  geometry  of  a  relationship,  consider  the  following  statement, 
water  is  pouring  cut  the  hole  in  a  bucket  faster  than  water  is  being  put  into  the 
bucket.'  Without  quantifying  either  flow  rate,  the  size  of  the  hole  or  the  capacity  of 
the  bucket,  we  still  have  a  firm  grasp  on  the  general  behavior  of  the  system.  As  Kane, 
et  al.  observed.  "Subjective  evaluations  generally  correlate  well  with  geometric 
understanding.  If  not  too  much  (precision)  is  asked  for  it  is  possible  to  get  more 
I  understanding  1."  [Ref  21;  p.  66]  Thus,  while  a  subjective  evaluation  is  relatively 
ambiguous  by  nature,  it  does  contain  useful  information  in  the  form  of  geometric 
relationships.  Further,  it  is  in  the  interpretation  and  evaluation  of  geometric 
relationships  that  mathematically  unsophisticated  people  can  use  their  intuition  and 
reason,  rather  than  rely  on  obscure  statistical  measures  and  overly  precise  numerical 
predictions. 

KSI.M  calculus  is  designed  to  im.part  a  feeling  for  linkages  that  cross  connect 
polmv  variables.  .As  a  prime  hvpothesis  we  assume  that  in  actual  policy 
implementation,  more  insight  is  needed  m  geometric  concepts  (the  connections  between 
■.  ariabics.  the  direction  of  ibrees.  and  the  tlmeshold  and  saturation  of  variables. )  Sucli 
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(.onsidcrations  liave  lar  more  ii!ipoirar.i.c  rh.m  arithmetic  speLil'ication  ol  paiamctcr'> 
Consequcatly.  KSIM  \s  dc'ogi'.cJ  to  impaii  an  appicciatioit  of  the  geometry  and 
dynamics  of  the  system  rather  tliait  an  appreciation  ot'  numerical  bookkeeping. 
|Ref.  21:  p.  Cj  With  these  notions  of  KSIM  philosophy  in  irand.  the  relative  softness 
sh  :nodcl  input  rcsiuircrnents  seems  to  make  more  sense.  I  he  required  inputs  are: 

•  rite  list  of  events  and  trends  considered  to  '•'c  tire  nuninuim  set  ersariabies  tb.at 
can  lulls  describe  the  svstom  under  studs, 

•  1  lie  initial  values  for  tltese  events  and  trends. 

•  riie  strength  and  mode  i.enltancing  or  inhibiting)  of  the  cross  impacts  between 
system  'variahles 

With  these  inputs  to  the  model  now  defined  and  the  output  of  the  model  with  its 
geometric  concept  described,  it  is  appropriate  now  to  develop  an  algorithm  to  relate 
tlic  two. 

C.  KSIM  MODEL  OF  AR.MS  TRANSF  ERS 

.As  a  review,  the  KSIM  model  has  the  following  five  properties: 

•  .All  system  variables  arc  bounded 

•  Variables  cliange  according  to  the  net  impact  of  all  other  variables 

•  S'ariable  response  to  net  impact  approaches  D  as  that  variable  approaches  its 
upper  or  lower  bound. 

•  ^'arlablcs  wall  c.xert  a  greater  impact  on  tlie  system  as  the  magnitude  of  those 
variables  grows  larger. 

•  ('omplcK  interactions  can  be  broken  down  into  binary  interactions. 

The  nratficmatica!  calculations  arc  carried  out  in  an  iterative  fashion.  With  At  being 
the  time  interval  for  one  step,  future  values  of  system  variables.  Xj(t  At),  arc 
computed  Irom  present  values,  X^fti.  according  to  equation  3.1 

rrdt) 

Xy.t  -  At)  =  Xj(t)  (eqn  3.1 ) 


1  he  term  7Tj(t)  is  derived  from  the  cross  impacts,  Y,jit). 

J 

JTj(t)  =  - 

1  -  AAt^bYijM'  *  Y,/t)i 


(cqn  3.2) 
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1  he  cross  impact  term.  Y|j(t).  is  a  Ilmctioti  of'  both  tiie  present  value  of  the  impacting 
variables  X.(t),  and  the  change  in  those  variables  dXj(t)'dt. 

YjjVt)  =  A-Xj(t)  -r-  B|-(dX-(t).dt)  (eqn  3,3) 

Kane  et  al  provide  an  e.xcellent  description  of  the  meaning  of  the  and  terms  in 
[Ref,  21j.  They  observe  that  the  A.^  term  describes  the  impact  variable  j  will  have  on 
variable  i  simply  because  of  its  e.xistcnce.  I'or  c.xample.  the  amount  of  sunshine  has 
this  type  of  impact  on  plant  growth.  On  the  other  hand,  the  elements  of  the  matrix 
describe  the  impact  that  a  change  in  the  value  of  variable  j  has  on  variable  i.  The 
impact  that  weather  changes  often  has  on  arthritis  pain  is  an  example  of  this  type  of 
impact.  ,A..  and  B..  mav  be  I'unctions  of  time  but  arc  almost  alwavs  constants.  (To 
require  the  user  to  estimate  a  functional  cross  impact  parameter  would  be  contrary^  to 
the  concept  of  simplicity  of  inputs’) 

In  Its  most  basic  form,  the  KSIM  model  is  described  by  these  three  equations. 
In  this  basic  form,  only  trend  variables  that  are  initially  constant  over  time  may  be 
used.  The  values  of  the  trend  variables  are  modified  according  to  the  cross  impacts  of 
other  trend  variables,  so  they  are  not  constant  in  the  presence  of  cross  impacts. 
Lipinski  and  Tydeman  sought  to  find  a  way  to  include  trend  variables  in  the  simulation 
that  were  not  initially  constant.  That  is.  the  trend  variables  of  interest  are  those  that 
take  on  new  values  over  time,  regardless  of  the  presence  of  cross  impacts,  fhese  initial 
nonconstant  \ariables  can  be  included  in  the  model  by  first  finding  the  function,  Q|(t). 
mch  that  it  describes  the  original  nonconstant  trend,  Xj(t),  in  a  recursive  way: 

Q,(t.  At) 

Xj(t  +  At)  =  Xdt)  (eqn  3.4) 

Now.  to  include  X(t)  in  the  KSIM  procedure  the  impacts  are  applied  according  to 
equ.aiion  3.5: 

Qj(t.  At)7r,(t| 

X.(t  At)  =  X.(t)  (eqn  3. .3) 


I 


s  s 


uc  now  note  that  any  given  (.uninlativc  distribution  (unction  (Cli)!')  can  be 


written  in  iterative  1'orni  as 


Q,(t.  ^tl 


C(t  -  Atl  =  C.(t) 


(eqn  3.6) 


'■\herc  Q^it)  IS  a  lunctioir  tliat  wiil  describe  the  CDF'  of  C^.  I’hus,  Q^tt)  is  a  function 
t'.iat  eields  the  a  priori  increment  to  the  CDF  at  each  time  step.  Our  objective  is  to  be 
able  to  learn  about  tire  beliacior  of  tlie  system  of  events  that  comprise  the  arms 
transler  process.  One  ol  the  dilliculties  of  using  the  KSI.M  model  is  to  define  e.xactlv 
u'hat  IS  meant  by  t.he  loosely  used  term  "system  variable”.  Therefore  we  shall  explicitlv 
'.leline  our  state  sariable  Satt,  the  value  ol  e\cnt  i  s  CDF  with  due  consideration  ol 
^ross  impacts,  by  the  following  recursive  relationship: 


Sat  ^  At)  =  S^(t) 


Q.(t,  At)7t_(t) 


(eqn  3.7) 


where  Q.(t)  is  the  a  priori  increment  to  event  i  s  CDF  at  time  t  and  7T.(t)  is  given  by 
equation  3.2.  lo  initialize  the  state  variable  Sj(t),  we  find  the  value  of  time,  (t^).  for 
which  C(ty)  =  0.0(.)l,  and  let  S(t^)  =  0.001.  In  the  absence  of  cross  impacts 

=  1),  equation  3.7  is  precisely  the  same  as  equation  3.6  and  the  system  variable 
lor  event  i  is  simply  the  CDF  for  i.  \Ve  shall  now  explore  the  properties  of  our  new 
model. 

D.  MODEL  PROPERTIES 

The  properties  of  this  new  model  are  examined  in  terms  of  Kane's  original  model. 
Kane's  first  property  states  that  system  variables  are  bounded  by  0  and  1.  Since  0.(t) 
generates  a  rnonotonically  increasing  function  (the  CDF),  it  must  be  less  than  1  over 
all  t.  We  are  given  that  Ttj(t)  is  greater  than  0  by  its  definition.  Since  S|(t^)  <  1  by  its 
definition,  then  the  value  of 


Similarlv,  since  S.ft  )  >  0, 


Thus  the  property  of  boundedness  in  the  open  interval  (0,1)  is  preserved  in  the  neu- 
model. 

Kane's  second  property  states  that  a  variable  increases  or  decreases  according  to 
whether  the  net  impacts  are  positive  or  negative.  When  there  is  a  net  positive  impact, 
rt.it)  <  1  from  equation  3.2.  Again,  we  know  that  Qj(t)  <  1,  so  the  product 
Q^it)7T-:tt  <  0.  Thus  we  can  positively  state  that  in  the  presence  of  a  positive  impact, 
S.'t)  will  increase.  However  if  the  net  impact  is  negative,  7t|(t)  >  1,  the  behavior  of 
S^i  t)  seems  to  depend  on  the  product  Q.(t)7t.(t).  Note  that  when 

Q.it)7t.(t)  <  1  S-(t)  increases 

^1-  r  r  ' 

Q^(t)7tj(t)  >  1  S.(t)  decreases 

Q.(t)7t;(tt  =  I  Sj(t)  remains  constant 

However,  what  seems  more  important  is  what  happens  to  the  state  variable  in  the 
presence  of  impacts  relative  to  what  happens  to  the  state  variable  without  impacts.  To 
explore  this,  consider  the  dilTerence  in  magnitude  between  S.(t)Qi^^'  and 
for  various  values  of  Tr.(t).  Observe  that  when: 

rt.(t)  <  1  Sj(t)Qi(^^  < 

K.ii)  >  1  s.(t)Qi^^^  > 

p  r  '  r  ^ 

From  the  above  relationships  it  can  be  seen  that  in  the  presence  of  net  positive  cross 
impacts,  the  state  variable  will  increase  above  its  a  priori  level  and  when  there  are  net 
negative  impacts,  the  state  variable  will  be  forced  below  its  a  priori  increment. 
Therefore,  we  can  conclude  that,  in  terms  of  our  model,  a  property  similar  to  Kane's 
second  property  is  satisfied. 

The  third  property  is  that  a  state  variable's  response  to  the  cross  impacts  of  other 
system  variables  will  decrease  to  0  as  that  state  variable  approaches  its  upper  or  lower 
bound.  To  investigate  this  property,  we  take  the  derivative  of  Sj  with  respect  to  time: 

S.(t  ^  At)-S.('t)  1] 

At  At 

I'rom  this  expression,  it  can  be  seen  that  as  Sdt)  approaches  0.  dSj(t)  dt  approaches 
/ere.  .Also,  when  S.(t)  approaches  I.  I|  approaches  0,  so  dS-(t)  dt 

SLii;  goes  to  U.  Clearly.  Kane  s  third  property  is  applicable  to  the  new  model. 
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Kane  s  fourth  property  is  that,  all  other  thing';  being  equal,  a  -.ariable  -.efl 
greater  inlluence  on  the  system  as  its  ntagnitude  increases,  hquaticn  s  '  st.it 
Y  ;it)  IS  a  function  of  both  the  niagnituue  of  S.it'.  and  the  rate  o!  lI. ar.ee  c 
fh.us.  as  the  magnitude  of  S^'t!  increases  it  ’.v;i!  produce  a  lareer  value  a;  ■/ 
conscsiuentiy  a  larger  net  impact  on  state  variable  i.  I  he  lourtli  KSIM  n^cperr. 

linallv.  ti'.e  I'lfth  assumption  is  that  onK  bin.tr.  interaL'icns  mav  '■^e  iv 
fit.s  property  holds  sm^e  svsrem  behavior  is  still  modelled  throuer  die  ^ross  , 


exert 

M  S  '  '  • 
t  >  and, 
holds 


eiven  bv  y  .t  ;i,  which  are  soieiv  a  function  of  event  pairs. 

^  '  .J 

I'hus.  our  new  model  for  the  value  of  an  event  C'DF'  with  considera:i,.'a'.  of  ^ross 
impacts  IS  trueiy  a  KSIM  type  cross  impact  procedure.  .Xs  such,  wc  are  n ''ouno 
tlte  necessity  to  find  a  suitable  numerical  .nterpretation  for  our  state  variable.  Sri.  our 
primap.'  concern  is  in  the  geometro  of  the  system  variables  and  the  elfects  the-,  have  on 
each  other  under  various  circumstances.  The  model  has  utility  regardless  of  cur  abilitv 


to  make  probability  statem.ents  !rom  the  values  of  its  state  variables  at  any  particular 
tinte.  since  inforntaticn  conveved  in  terms  of  probabilities  is  much  less  appealing  than  a 
visual  display  of  the  growth  rates  and  magnitudes  of  the  state  variable  over  time. 
.-\gain  we  refer  to  the  basic  philosophy  of  the  KSIM  procedure:  The  significant 
diference  in  orientation  between  our  procedures  and  most  other  methods  is  that  we 
emphasize  the  geometrx'  of  the  system,  the  structural  relationships  between  the 
variables  wiule  standard  procedures  tend  to  emphasize  arithmetic  details,  the  precise 
'pecification  of  coelFicients  and  parameters."  (Ref.  22:  p.  2S6] 

Now  we  si'.all  take  up  the  task  of  implementing  this  model.  The  first  part  of  this 
umk  ..oriMsts  of  two  elements;  the  development  of  an  expression  for  Qat)  and  providing 
'c-me  meth.od  to  estimate  the  parameters  needed  to  specify  QM).  The  rest  of  this 
vf...p'.er  adJrevves  the  details  of  this  first  task. 


E.  DETERMINATION  OF  QCH 

I  he  arms  transfer  process  can  be  broken  down  into  eight  discrete  events.  The 
quev'.ion  is  what  type  of  CDF'  should  we  use  to  model  these  events.'  Mitchell  Bloom 
:r..;ke'  a  strong  case  for  using  the  logi'itws  equation  to  describe  the  CE)F'  of  an  event. 
:Rcf  12:  p.  IS5] 

Ihe  probabi'itv  density  function  'PDFb  for  the  logistics  distribution  ix  an  inverted 
I  -xhape  Mhgure  .'.lx  The  interpretatmn  of  this  PDF  is  that  there  is  a  negligible 
pr:b.;'v;,rv  ,.'f  i.xLurrence  .u  e.irlv  timex  vvhwii  rises  to  a  maximum  and  fails  back  olf  :o 
1  :.cc'..g;:':e  .eve!.  I  lie  concept  of  representing  event  probability  denxitiex  ;n  tins 


irannor  is  neither  strange  nor  new.  It  is  oltcn  saiJ  that  il  an  event  Joes  not  occur 
belorc  a  certain  point  in  time,  its  probability  Jeciincs  rapidly  until  it  is  almost  certain 
not  to  occur.  Lxamples  include:  waiting  for  a  rricnd  to  arme  at  a  predetermined  place 
and  tune  or  waiting  lor  the  partv  one  has  called  to  answer  the  telephone.  IRcf'.  12:  p. 
IS'j  Ihis  notion  is  fully  compatible  with  describing  arms  transfer  events,  an  order  for  a 
o  capon  scstcm  has  a  ma.ximum  probabilitc  of  occurrence  some  days  alter  a  negotiation 
.mJ  il’th.e  order  is  not  placed  within  a  certain  period  of  time,  it  probably  will  not  occur, 
l  or  tliese  and  for  reasons  of  mathematical  tractabihty  m  the  K.SIM  procedure,  the 
losistics  curve'  will  be  used  to  model  arms  transler  esent  ('DF  's. 
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Il  It  IS  liter  deckled  tli.at  s.'me  ''ihcr  distribution  is  more  appu'priite 
ue  used  to  lind  C),'t)  m,iv  be  used  .  n  .ii:-.  diilerenlabic  liinction. 


1  Event  Parameter  Equations 

The  logistics  CDF  is  given  by  equation  \.logeq..  (For  ease  of  notation  anJ 
u.i  Jerstanding,  the  remainder  of  this  development  -.vill  consider  only  one  event  and  will 
omit  its  subscript,  i.  When  a  final  computational  form  is  reached,  the  subscript  will  be 
replaced.  i 


C  = 


1 

1  -  e\p(  -  at  -  P) 


<  eqn  3.9  i 


'.vhere  C  is  the  value  of  the  CDF  at  time  t,  and  a  and  P  are  constants  to  be  determined. 

To  compute  a  and  p,  ue  estimate  the  time  at  which  the  cumulative  probability 
IS  equal  to  O.b  denoted  t^.  The  time  at  which  the  cumulative  probability  is  equal  to  o.l 
will  be  denoted  t^.  iSee  Figure  3.2)  Thus, 

'.vhen  ^  “  hi  •  ^  ~  “ 

t  =  t^  ,  c  =  0.1 

Equation  3.9  can  be  written  in  logarithmic  form: 

-at  p  =  !n((  1  C)  -  1)  (eqn 


Substituting  the  values  of  C  for  t^  and  t^  shown  above  into  equation  iSilnlogcdf  and 
scl'.'ing  for  a  and  P  yields; 


a  =  (I  (tjj-t^i!ln9  (equ  3.11) 


P  =  't^  ( t|^  -  t^ Mln9  ieqn3.12' 


order  to  obtain  the  fuiKtion  Qai  Irom  the  iterative  equation  3,^,  we 


.itr-  e 


of  (  u  ith  re'ipect  to  t; 


Since  we  want  to  get  out  ot  the  recursive  equation  and  find  a  closed  form  solution  for 
Qit).  we  will  set  the  right  hand  side  of  equation  3.15  equal  to  Cit)^''’'*: 


C(t)Q-^'  =  C(t) 


a(  At  lexpi  -  at  -  P) 


1  -  expi  -  at  -  p ) 


] 


(eqn  3.16) 


1  ina'.lv.  taking  the  natural  logarithms  of  equation  3.16  yields  equation  3.1"', 


0(t)InC(  ti  =  InCit)  in 


[- 


a(At)exp(  -at  ^  p) 


1  -  expi  —  at  -  pi 


(eqn  3. r 


By  diMU.ng  equation  3.1'  by  InC(t),  wc  have  an  expression  for  Q(t).  The  subscripts  are 
now  returned  and  our  final  closed  form  equation  for  Q.(t)  is 


In 


Q.it)  =  1  - 


''  a.(  J 

1 - ^ 

1  - 


.( Atie.xpi  —  at  +  p. 


expi  -  au  -f-  P^ ) 


r] 


In 


[_  1  -  expi  -a.t  -  p.) 


(eqn  3. IS) 


F.  E^  E.NT  P.ARAMETER  ESTIMATION 

In  order  to  implement  this  model  the  user  will  be  required  to  supply  estimates  for 
the  lollowmg  parameters  for  each  event: 

•  Time  when  CDF  is  equal  to  0.5.  t^ 

•  Time  when  CDF  is  equal  to  0.1.  t^ 

•  Strength  and  mode  of  the  existance  of  this  variable  on  every  other  variable 

•  Strength  and  mode  of  an  increase  in  this  variable  on  every  other  variable. 

If  the  user  is  very  familiar  with  the  arms  transfer  process,  estimation  of  these  values 
mav  net  prove  to  be  difficult.  .As  John  Mather  points  out  though.  ".A  more  useful 
de\elopment  would  be  the  provision  of  ■^ome  means  of  estimating  the  initial  \alues  of 
th.e  'e\ents)  from  available  data.'  [Ref  23:  p.  21]  If  a  first  estimate  of  event  parameters 
..ould  be  rroMded  by  an  existing  data  base,  then  this  would  greatly  reduce  the  burden 
:n  tl’.c  anal;.  St  and  proxide  a  solid  basis  for  discussions  regarding  the  acceptance  or 
alteration  of  the  estim.ated  parameters.  Third  Point  Systems.  Inc..  Monterey. 
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California  has  made  a  substantial  effort  to  record  arms  transfer  data  for  almost  ever>- 
countin'  in  the  world.  It  is  possible  to  extract  event  parameter  estimates  (with  the 
exception  of  the  cross  impact';i  from  this  data  set. 

1.  Event  Data  Set 

The  .Arms  Transfer  Data  Set  was  developed  by  Third  Point  Systems  to  record 
individual  arms  transt'or  events  which  could  be  aggregated  to  describe  the  entire  arms 
transfer  process.  The  purpose  for  the  creation  of  the  data  set  was  to  aid  foreign 
policy  decision  makers  to  evaluate  the  patterns,  purposes  and  effects  of  international 
arms  transfers."  iRef.  24;  p.  5]  The  other  major  sources  of  arms  transfer  data,  such  as 
SIPRI  and  .ACD.A.  have  concentrated  only  on  the  value  of  the  weapons  translerred, 
and  the  date  the  transrer  took  place.  .Adding  further  to  the  confusion,  these  sources 
use  diiTerent  dates  for  the  actual  arms  transfer;  some  use  the  date  the  contract  is 
signed,  seme  use  tiie  date  the  order  is  placed,  and  others  use  the  actual  date  of  deliven.'. 
However,  an  arms  transfer  consists  of  many  discrete  events  and  focusing  only  on  the 
cash  value  of  a  weapon  and  its  date  of  transfer  severly  limits  a  thorough  analysis  of  the 
arms  transfer  process.  Recognizing  this  unneccesary  limitation.  Third  Point  Systems 
began  to  build  an  extensive  data  set  that  recorded  each  arms  transfer  event.  With  this 
data,  the  analyst  could  interpret  the  subtle  changes  in  the  attitudes  of  various  nations 
towards  '.veapons  sales  and  procurement.  Having  done  this,  the  foreign  policy  analyst 
nught  be  better  prepared  to  attempt  to  influence  other  country's  actions  and  provide  a 
small  measure  of  control  over  the  arms  transfer  process. 

Third  Point  Systems  identified  fifteen  distinct  events  that  comprised  the 
various  stages  of  an  arms  transfer.  For  the  purposes  of  this  study,  some  of  these 
events,  such  as  capture,  were  eliminated  as  not  being  particularly  influential  on  the 
-I'ms  transfer  process.  Other  events,  such  as  reject,  refuse  (recipient)  and  reject  refuse 
''Upphen.  were  combined  into  one  event  since  they  were  not  different  enough  to  be 
considered  separately.  Thus,  the  event  set  used  in  this  analysis  consists  of  the 
fcilowing  eight  events  (the  descriptions  are  based  on  the  .Arms  Transfer  Handbook 
jRef  24):j 

•  Meet/Visit.  This  event  occurs  when  two  countries  meet  to  consider  an  arms 
transfer  issue.  It  includes  the  exciiange  of  notes,  messages,  and  information  as 
well  as  the  actual  meeting  between  members  of  the  two  countries. 

•  Propose/ Request.  When  a  supplier  country  offers  a  weapons  system  or  support 
package,  or  when  the  recipient  directs  a  request  to  the  supplier,  then  the  event 
IS  coded  as  a  Propose  Request. 
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•  Evaluate/Negotiate.  Events  in  this  categon-  indicate  that  a  countrv'  is  in  ’he 
process  of  considering  the  purchase  of  some  svitem.  Evaluations  can  be  both  a 
supplier  countrv'  considering  a  sale  or  a  recipient  countrv'  considering  a 
purchase.  Also,  agreements  between  industries  or  governments  to  produce  a 
specific  piece  of  equipment  for  development  is  considered  an  Evaluation, 

•  Reject/ Refuse.  This  event  occurs  when  either  the  recipient  or  supplier  reacts 
negatiN'ely  and  rejects  its  counterpart's  otTer  or  proposal. 

•  Order.  If  a  contract  is  signed  or  awarded,  or  if  an  agreem.ent  is  made  to 
purchase  or  coproduce,  then  the  event  is  coded  as  an  Order. 

•  Delivery.  This  type  of  event  includes  ^i.e  actual  delivery  of  a  weapons  system, 
tile  return  to  an  operational  state  after  overhaul  by  an  e.xternal  supplier,  and 
the  licensed  production  or  coproduction  of  weapons  systems. 

•  Increase.  Decisions  to  resume  a  previously  halted  deliver}',  reduce  sanctions  or 
supplement  a  previous  order  are  all  considered  Increases. 

•  WitlidraNN/Cancel.  Events  in  this  category  are  those  where  the  supplier  or 
recipient  slows  production,  reduces  the  quality  of  weapons,  or  cancels  an 
agreemient. 

Eigure  3..'  shows  the  relationship  between  these  events  and  possible  linkages  that  e.Kist 
between  them. 

2.  Data  Set  Coding 

The  ev  ent  data  record  is  made  up  of  two  parts:  the  first  is  an  analytical  section 
with  a  strict  coding  scheme,  and  the  second  is  a  narrative  summary  of  the  event.  The 
iirst  section  contains  the  data  which  can  be  used  to  provide  estimates  of  event 
parameters.  This  section  will  be  discussed  in  detail,  .-k  sample  coding  of  an  arms 
transfer  event  appears  in  Table  1. 

The  analytical  section  of  the  data  will  now  be  broken  down  line  by  line. 

Line  1  Field  1:  .Month  ( 1 ) 

Field  2:  Day  ( 31 ) 

Field  3:  h'ear  (86) 

The  information  contained  on  line  1  is  self  explanatory.  Event  dates  are  only 
i^oded  ’.vhen  the  event  occurs;  if  an  event  is  predicted,  it  is  coded  as  a  Comment. 
Sirralarly,  if  an  event  is  alluded  to  which  occurred  more  than  six  months  prior  to  the 
source  date,  it  is  not  coded;  it  is  assumed  that  these  events  were  coded  earlier. 

Line  2  Field  1:  .Actor  country  (.A LG  =  .Algeria) 

Field  2;  Domestic  actor  (GVT  =  government) 

Field  3:  Event  code  (53  =  Request  Inquire) 

Field  4:  Domestic  target  iGVT  =  government) 


REJECT 

REFUSE 


ORDER 


TABLE  1 

EVENT  DAEV  RECORD 


31  66 


53 


GVT  BRA 


EE-9 


AVREWH  El 


Algeria  shows  interest  in  the  Engesa  Cascavel 
(Brazil)  EE-9  armored  combat  vehicle  with  a 
9Cm.n  g'cn. 


Field  5;  Target  country  (  BR,\  =  Brazil) 

An  actor  country'  is  the  countin'  who  initiates  an  arms  transfer  event.  The 
actor  countrv'  can  he  either  the  supplier  or  the  recipient;  the  key  being  which  countin' 
initiated  the  action.  For  example,  in  a  Request,  tlie  recipient  is  the  actor,  whereas  in 
an  Oi'er,  the  supplier  is  the  actor.  The  target  countin'  i^  the  object  of  the  action 
initiated  b;.  the  actcr  ccunirc.  The  domestic  actor  and  the  domestic  target  refer  to  the 
pervcn.  organization  or  group  within  the  respective  ccuniin  tl'.ai  i'^  respon‘;’bie  I'or  the 
event.  If  n.o  domestic  actor  target  is  specified,  the  government  a^'mmed.  Fb.e  event 
..ode  refer':  to  the  fi.Teen  specific  events  that  comprise  the  arn'is  tran>;er  proce-:';. 

Line  3  Field  1:  Terms  of  sale  (3.''  =  other) 

The  terms  of  sale  of  a  weapons  svstern  is  indicated  if  the  inform.ation  is 
available.  Examples  of  these  terms  include  coproduction,  clNets,  gifts  and  credit. 

Lined  Field  1:  3^'eapon  s\ '-tern  >  EE-9 ) 

Field  2:  Equipment  typciAVREWH  El) 

ihe  weapon  system  is  the  broad  category'  of  hardware  support.  The 
j^iaipment  trne  uses  a  code  to  snccificaiiv  break  down  the  svstem  to  the  component 


AUer  the  data  is  collected,  it  is  assembled  in  the  form  of  "stonhnes'  for 
review.  The  storvTine  simply  collects  all  events  for  a  given  countrv-  both  when  it  is  an 
actor  and  when  it  is  a  target.  In  this  way,  the  reviewer  can  see  if  there  exists  a 
coherent,  logical  progression  of  events.  Once  the  reviewer  is  satisfied  wuh  the  recorded 
data,  It  IS  loaded  into  a  mass  storage  device. 

'  Parameter  Estimation  Technique 

To  obtain  an  estimate  for  t^  ithe  time  when  the  CDF  is  equal  to  0.5)  and  t^ 
■the  time  when  the  CDF  is  equal  to  0. 1),  one  must  find  a  series  of  events  between  two 
countries  that  deals  with  the  same  weapon  system.  For  example,  the  olFer  of  the 
L  nited  States  to  seil  F  .A-IS  s  to  Israel,  the  evaluation  by  Israel  and  the  subsequent 
order,  etc.  rraght  be  used  as  a  source  of  one  piece  of  data.  The  number  of  days 
between  event  pairs  is  recorded  for  all  storylines  that  exist  in  the  data  set  for  similar 
weapons.  (Presumably  it  takes  longer  to  evaluate  the  purchase  of  a  squadron  of 
I  .A-FST  than  a  load  of  M-16  rifles.)  A  CDF  would  then  be  constructed  from  this 
data.  The  day  when  the  CDF  was  equal  to  0.5  would  be  tj^  and  the  day  when  the  CDF 
was  equal  to  0.1  would  be  t^.  This  approach  requires  a  great  deal  of  data  since  the 
number  of  storx'iines  in  the  data  set  is  relatively  small,  however,  this  analyst  only  had 
access  to  a  small  portion  of  the  data  set.  so  an  extensive  study  in  this  field  was  not 
possible.  (Indeed,  this  type  of  study  would  be  expansive  enough  to  support  an  entirely 
separate  research  report!)  In  the  data  provided,  seven  samples  were  found  that 
contained  the  linkage  between  Evaluation  Negotiation  and  Order  on  similar  pieces  of 
equipment.  The  parameter  estimation  technique  will  be  demonstrated  with  these  seven 
d.ita  points.  The  raw  data  are  shown  in  Table  2.  These  seven  data  elements  were 
found  in  storylines  from  various  countries.  Obviously,  to  implement  this  technique, 
one  would  restrict  themselves  to  the  arms  transfer  data  from  the  countin'  of  interest. 

To  construct  the  CDF,  we  compute  the  probability  that  the  number  of  days 
between  an  evaluation  and  order,  a  random  variable  T,  will  be  less  than  or  equal  to  the 
various  elapsed  times,  t,  from  Table  1  That  is. 

Pr(T  <  t) 

The  results  of  this  calculation  arc  shown  below. 


t 

60 

64 
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DATA  FOR  FSlIMATkjN  OF  FVFNT  PAR.AMFTFRS 


of 


iuation  Negotiation 

Date  of  Order 

ITapsed 

Ii9  ii2  S5 

11  16  85 

> 

IIS  1)1  85 

no  3()  s5 

6() 

in  28  85 

12  31  85 

65 

"2  ns  86 

u3  o6  86 

in  n5  85 

12  10  85 

66 

01  2n  86 

04  19  86 

S9 

')'  01  85 

(,)2  r  86 

•>  *; 

Theroi'cre.  '.ve  '.voula  use  66  as  our  estimate  oft,  and  2’’  as  our  estimate  oft  .  .At  this 

n  e 

point,  we  rccogni/e  that  with  all  eight  e\ent  CDF's  going  from  <)  to  1.  the  plot  will  get 
eiuttered  very  quicklv.  Further,  we  know  that  by  multiplying  the  CDF  by  a  con'^tant 
Joes  not  ruin  any  of  the  five  properties  of  the  model  previously  discussed.  We  could 
arbitrarily  select  values  for  these  constants  since  they  have  no  numerical  m.eanmg  in 
our  ntcdei.  but  to  give  one  mere  piece  of  visual  data,  we  will  use  the  relative  frequency 
of  each,  event  in  the  Third  Point  data  svstem  and  denote  this  constant  M. 


I  ^  exp.  -  ut  *  P) 


(eqn  3. FT 


,\s  an  example,  suppose  Saudi  .Arabia  was  the  countn.’  in  question.  Table  3 
■'ummarizes  Saudi  .Arabia's  arms  transfer  activity  for  the  year  1985  -1986.  From  this 
'.;mm.;r’.  it  c^r.  be  '•ecn  that  the  relatise  frequency  of  Order  is  0.4"4.  so  this  would  be 
tne  -.aiue  of  M. 


TABLE  ? 

SAUDI  AR.ABIA  ARMS  TR-\NSFER  ACTIVITY  ( 1985  -  1986) 


Event 

Nr  ofOccurances 

Percentage 

Propose  Offer 

16 

5.8 

Reuuest 

r 

6.2 

Evaluation  Negotiation 

■() 

25.5 

Reject  Refuse  (supplier) 

I.l 

Reject  Refuse  i recipient) 

3 

I.l 

Order 

130 

47.4 

Delis  er 

13 

4.7 

Increase 

17 

6.2 

\5'ithdraw  Cancel 

"> 

0.7 

As  stated  earlier,  the  data  necessar\-  to  make  reasonable  estimates  for  all  event 
parameters  was  unavailable.  The  values  of  the  parameters  used  in  further  analysis  are 
based  on  the  author  s  estimates  and  in  no  way  are  intended  to  represent  any  one 
ccuntiT  or  arms  transfer. 

G.  KSIM  ALGORITHM 

To  be  effective  as  a  policy  analysis  tool,  the  KSIM  model  for  arms  transfers  must 
be  programmed;  the  calculations,  while  not  dilTicult,  are  too  numerous  to  be  handled 
manually.  .Additionally,  the  KSI.M  program  should  be  portable  so  that  it  is  not  bound 
to  a  main  frame  facility  and  become  relatively  inaccessible  to  the  intended  user.  Tor 
these  reasons,  it  was  decided  to  program  the  KSIM  model  on  an  IBM  XT  personal 
cemputer.  The  program  follows  the  algorithm  outlined  below: 

Step  1.  Input  the  event  parameter  estimates  t^.  t^  and  M  for  all  events. 

Step  2.  Input  tiie  cross  impact  matrices. 


1%' 


I 

i'  • 


Step  3.  Compute  the  values  of  a  and  P  according  to  equations  3.11  and  3.12. 

Step  4.  Compute  the  time  when  the  event  cumulative  probability  is  equal  to  O.nOl.  If 
this  time  IS  less  than  0,  then  compute  the  value  of  the  cumulative  probability  at  time  0. 

'  L  '^e  equation  3. 19 1. 

Step  5.  Deternune  the  start  time  of  the  simulation.  The  start  time  will  he  the 
minimum  of  e,wh  e\ent  s  t^  or.  if  this  is  less  than  0.  the  start  time  will  be  0.  Set  the 
simulation  time  equal  to  the  start  time. 

Step  6.  Check  all  events  to  see  if  their  t^^  is  less  than  or  equal  to  the  current  simulation 
time.  .Any  such  events  become  part  of  the  simulation  with  an  initial  value  of  O.t'Ol. 

Step  7.  Compute  y  ■  according  to  equation  3.3. 

Step  8.  Compute  n.  according  to  equation  3.2. 

Step  9.  Computes, It  -  At)  according  to  equation  3.7. 

Step  10.  Repeat  Step  6  until  the  stop  time  (selected  by  the  user)  is  reached. 

This  algorithm  was  programmed  using  Turbo  Pascal  (version  3.01a).  .A  listing  of  the 
source  code  appears  in  the  .Appendi.x.  To  make  the  program  functional  as  an  analysis 
tool,  it  had  to  be  user  friendly,  interactive  and  allow  various  parameters  to  be  changed. 
.As  such,  the  program  is  menu  driven,  and  allows  data  editing  of  event  parameters  and 
cross  impact  values.  Upon  termination  of  the  session,  the  current  values  of  all  event 
parameters  and  the  cross  impact  matrix  are  stored  on  the  disk.  When  the  program  is 
loaded,  It  recalls  these  parameters  from  the  disk  and  returns  the  system  to  its  state  at 
tb.e  end  of  the  previous  session.  The  next  chapter  will  show  some  results  from  this 
program  and  discuss  how  the  model  might  be  used  to  investigate  various  arms  transfer 
policy  decisions. 


tv: 


IV.  RESULTS  .AND  CONCLUSIONS 

A.  MODEL  RESULTS 
i.  Input  Parameters 

.As  stated  in  the  previous  section,  the  data  used  in  the  testing  of  the  KSIM 
arms  transier  model  were  estimated  by  the  author.  Table  4  shows  the  initial  event 
parameters  used  in  the  simulation.  The  values  tor  the  relative  frequency  of  occurrence. 
.\1,  came  from  Table  3.  Since  there  was  no  value  for  the  event  .Meet  Visit,  U.ul  was 
a^ed.  The  two  e\'ents.  Propose  Offer  and  Request,  were  lumped  together  to  Ibrm  one 
e\ent,  Propose  Request  (see  the  subsection  on  Event  Data  Set).  The  value  of  \I  for 
Propose  Request,  then,  is  the  sum  of  the  values  for  Propose  Offer  (0.058)  and  Request 
■  o.'  icd  j. 


TABLE  4 

INPUT  PAR.AMETERS  FOR  KSl.M  SI.MULATION 


Event 


Relative 

Frequency 


.Meet  \'isit 

20 

10 

.0100 

Propose  Request 

35 

15 

.1201 

Evaluate  Negotiate 

55 

40 

.2554 

Reject  Refuse 

60 

50 

.0220 

Order 

80 

6() 

.4  •’44 

Deliver 

100 

85 

.04  "4 

Increase 

130 

100 

,0602 

Withdraw  Cancel 

140 

120 

.00'’0 

Ihe  i’.rst  run  ot  the  nioJel  used  a  value  of  0  for  all  A.,  and  B...  This  implies 

;j  i; 

that  'here  are  no  cress  impacts  and  the  results  should  yield  a  truncated  (1DI‘:  tf.e 
salue  ct'the  CDF  multiplied  b>  its  relatise  frequenc;'. .  As  can  be  seen  iVcm  f  igures  4  1 
.m.d  -  2.  this  IS  indeed  the  case. 

rite  cross  impact  matrices  were  then  estim.ated  using  Figure  3.?.  If  there  uas 
„  i.iie  .  i.ite^t.iig  f.'.  o  events,  it  wa^  assum.ed  that  an  increase  in  the  probabilitv  of 
o.aurreit^e  tite  f.rst  will  have  seme  etfect  on  the  probability  of  the  second.  Recall 
tit. It  positive  '..hues  are  enhaitwing  impacts,  while  negative  values  are  inhibiting.  The 
'...lie  .wed  for  tite  derivative  i  impacts  was; 

0  =  None 

1  =  Mild 

2  =  Strong 

3  =  Overwhelming 

fai'ie  '  sitows  the  values  for  the  derivative  cross  impact  matrix  that  were  used  in  the 
s;ntt.i..t.cit. 

I  ite  d-ta  for  tite  constant  cross  impact  matrix  (  A|.)  was  estimated  in  a  simalar 
nt.inner.  Recall  that  this  matrix  is  the  impact  that  the  first  variable  has  on  the  second 
'lit. pi’.  be..ause  of  its  existence.  For  this  matrix,  however,  the  following  scale  was  used: 

0.0  =  None 

(1.1  =  .Mild 

0.2  =  Strong 

(1.3  =  Overwhelming 

fa'V.o  ('  displays  the  values  used  for  the  constant  cross  impact  rriatri.x: 

rite  data  were  entered  into  the  model  and  a  period  of  2(.30  dav  s  was  simulated. 
Tite  rc'ults  of  this  first  run  appear  in  Figure  4.3.  The  output  shows  that  the  cross 
impacts  have  a  great  influence  on  the  events  Evaluate  Negotiate  .  Reject  Refuse  and 
Order.  These  three  event  state  variables  reached  their  upper  bounds  even  though  their 
lii.al  values  without  cross  impacts  was  very  low.  To  understand  vvhv.  we  must  closely 
examine  the  cross  impact  matrices  and  couple  the  values  we  see  in  the  tables  vvitli  the 
i'ei'.avmr  we  observe  on  the  plot.  .\s  an  example,  let  us  try  to  determine  whv  the  state 
v  ari.f''ie  nor  Reye.;:  Refuse  grew  so  rapidly  and  to  such  a  large  magnitude. 


TABLn  5 

DtRIV ATIVL  CROSS  IMPACT  MATRIX 


IMPACTING  EVENT 


B 

:j 

1 

> 

4 

N 

6 

8 

1 

i) 

0 

0 

0 

() 

0 

0 

0 

1 

0 

0 

0 

0 

0 

0 

(.) 

() 

3 

0 

i) 

0 

0 

0 

0 

I.MPACTED 

4 

0 

0 

0 

.2 

0 

0 

0 

EVENTS 

N 

0 

0 

1 

-3 

0 

0 

0 

0 

6 

0 

0 

0 

0 

-y 

0 

0 

0 

- 

I) 

0 

0 

0 

1 

1 

0 

_*> 

8 

() 

0 

0 

0 

1 

1 

-2 

0 

E\cnt  I:  Meet  Visit 
E\cnt  2:  Propose  Request 
Event  3:  Evaluate  Negotiate 
Event  4:  Reject  Refuse 


Event  5;  Order 
Event  6;  Deliver 
Event  7:  Increase 
Event  8:  Withdraw  Cancel 


Refer  to  Tables  5  and  6.  Reject  Refuse  has  a  positive  derivative  impact  from 
Evaluate  Negotiate  and  a  negative  impact  from  Order.  Thus,  when  there  is  a  positi\e 
slope  in  Evaluate.  Negotiate,  the  incremental  increase  in  the  state  variable  will  enhance 
the  growth  of  Reject  Refuse.  Similarly,  as  the  state  variable  for  Order  increases,  it  will 
inhibit  the  sro'Ath  of  Reject  Refuse.  Now  since  the  B..  values  are  the  same  (2)  the  net 
derivatr.  e  impact  will  depend  on  the  dilTerence  in  slope  between  the  two  impacting 
sariablcs.  It  ^an  also  be  seen  that  the  state  variables  for  Propose  Request  and 


TABLn  o 

CONSTANT  CROSS  IMPACT  MATRIX 


I  M  P  A  C  T 1  N  G  EVENT 


A.. 

'•1 

1 

-» 

A 

4 

N 

6 

s 

1 

1) 

t) 

(1 

0 

0 

1) 

0 

0 

-V 

-1 

n 

(1 

() 

1) 

IJ 

(t 

; 

i) 

•> 

(> 

(1 

() 

0 

0 

0 

IMP.XCTED 

4 

(  ) 

.1 

0 

(,■) 

0 

0 

0 

EVENTS 

N 

() 

.1 

-3 

0 

0 

(J 

0 

6 

o 

«) 

1) 

,  2 

A 

0 

0 

1) 

- 

(■) 

0 

0 

*  • 

.1 

0 

-.1 

8 

(.1 

0 

0 

.1 

.1 

.1 

-.1 

0 

Event  1:  Meet  Visit 
Event  2:  Propose  Request 
Event  C  Evaluate  Negotiate 
Event  4;  Reject  Refuse 


Event  5:  Order 
Event  6:  Deliver 
Event  7:  Increase 
Event  S:  Withdraw  Cancel 


Evaluate  Negotiate  have  positive  constant  impacts  on  Reject  Refuse.  Whenever 
E'. aluate  Negotiate  and  Propose  Reque'-t  have  any  magnitude,  they  will  enhance  the 
growth  of  Reject  Refuse.  Now  look  at  Figure  4.1.  We  can  see  th.at  both 
Propose  Request  and  Evaluate  Negotiate  have  rapidly  increasing  state  \anables  (large 
slope -and  tliat  their  grc.nh  is  both  faster  (higher  slope)  and  of  greater  magnitude  than 
Order  t.p  to  da\  number  During  the  period  from  day  0  to  "5  then.  Reject  Refuse 
i.as  net  rositr. e  impact  that  accelerates  its  crowth.  In  fact,  the  erowth  rate  is  so  hieh 


v: :;'.o  KSIM  prc^eJare. 

Lo:  us  no'A'  Lor.siJcr  the  elTcct  c!' a  large  positive  impact  on  a  state  variable. 
!:  m.i;.  '•'c  cur  behef  that  several  events  have  a  positive  constant  cross  impact  on 
Del. '.cr,  so  we  shall  ^et  .Aj^  i  Meet  Visit  on  Deliver)  equal  to  0.1, 

■ alu.ite  Negotiate i  equal  to  h.2.  and  change  .Ajj^  (Reject  on  Deliver)  from  -0.2  to  0. 
Ir.e  model  is  run  for  a  200  day  simulation  and  the  results  appear  in  Figure  4.5.  .Again 
■■ve  see  th.it  the  net  positive  impacts  did  have  a  marked  elTect  on  the  state  variable 
Deliver,  but  the  secondatp.'  elTect  of  greatly  enhanced  grouih  of  Increase  is  even  mere 
;mor..' tinned,  T'ne  caused  of  this  secondary  effect  can  be  traced  to  the  entn.'  B^^-i  the 
.r.wre.l^;:■.g  value  cf  Deliver  has  a  positive  derivative  impact  on  the  state  variable 
In^reave.  \5'e  continue  to  modify  the  cross  impact  parameters,  observing  the  elfect  of 
these  Lhanges  until  we  obtain  a  set  of  entries  in  both  of  the  cross  impact  matrices  that 
result  in  '•vstem.  behavior  that  appears  correct.  It  is  during  this  process  of  refining  the 
nmdel  that  'he  analyst  will  learn  a  great  deal  about  the  interplay  of  state  variables 
.v.:h:r.  the  system.  .Many  users  of  the  KSIM  procedure  state  that  this  acquired 
knc'.vledge  is  often  of  more  value  than  the  results  of  various  simulation  runs.  It  will  be 


.;^"umed  that  the  state  of  the  model  as  'specified  by  the  cross  impact  matrices  in  Tables 
md  0  is  satisfactory  at  this  point,  and  proceed  with  investigating  some  of  the  arms 
■rui'ler  policies  suggested  in  chapter  one. 


2.  Promoting  the  Transfer  of  Defensive  Weapons 

Chapter  oite  al’uded  to  the  fact  that  some  weapons  lead  to  greater  political 
stability  than  others  by  their  technical  nature  and  defensive  mission,  and  that  it  m.a\  be 
desirable  to  promote  the  transfer  of  these  weapons.  One  policy  to  achieve  this  goal 
might  be  to  actively  "sell'  the  weapon  to  target  countries  by  demonstrating  its 
capabilities  and  proposing  its  purchase.  We  can  now  check  to  see  if  this  policy  will 
ha'.e  tile  desired  etfect.  To  do  this,  the  value  of  M.  the  relative  frequency  of 
occ..rrence.  will  be  changed  from  0.1201  to  O-.'iOO  for  the  event  Propose  Request.  The 
meaning  in  terms  of  probabilities  or  CDF's  of  increasing  .VI  from  0.1201  to  O.?0  is  not 
dear,  but  we  only  wish  to  model  an  "increased  effort"  in  this  state  variable  and 
in.-rcasing  the  value  of  .M  tends  to  convey  this  idea  of  "increased  effort".  The  output 
from  this  20f)  day  simulation  are  shown  in  Figure  4.6 

The  results  of  this  run  dem.onstrate  that  there  is  no  perceptible  change  in  any 
of  the  system  variables  due  to  the  increased  elTori.  In  fact,  observe  that  even  the  state 
variable  Propose  Request  is  unalTected.  The  value  of  M  was  subsequently  increased 
several  times,  but  the  results  were  always  similar  to  that  displayed  on  Figure  4.6.  We 
might  conclude  that  the  interactions  between  variables  themselves  are  more  important 
in  the  outcome  of  the  simulation  than  are  the  mitial  starting  points.  Can  we  draw  any 
reiesai'.t  policy  decision  information  from  this  simulation.^  Certainly  one  would  not 
state,  based  on  the  results  of  this  run.  that  an  increased  effort  in  Proposals  has  no 
beneiit.  Hosvever.  the  model  does  tend  to  indicate  that  the  arms  transfer  process  is 
quite  resistant  to  change  from  simple  event  parameter  changes. 

3  E.xtended  Congressional  Debate 

.Another  arms  transfer  policy  proposed  in  chapter  one  is  that  the  arms  transfer 
process  should  be  made  more  "viscous"  by  extending  the  length  for  Congressional 
debate.  \\'hat  might  the  effect  on  arms  transfers  be  if  this  policy  were  implemented? 
fhe  parameters  for  Evaluate  Negotiate  are  currently  set  at 

t^  (time  when  CDF  =  0.5)  =  55 
tg  (tim.e  when  CDF  =  0.1)  =  40 

To  simulate  extended  Congressional  debate  these  values  can  be  modified  to 
tj^  (timie  when  CDF  =  0,5)  =  75 


DAYS 


Again  the  ntoJel  >  ru: 
e\ten^-»eL»  0 ong^ essiona; 
probes'.  In  faet.  ti.e  a: 
,0''  'teen  is  to  oeia'.  t; 


-r  2""  Ja'.  s.  3y  examining  Figure  4."'  it  can  be  seen  that 
'ate  u'es  not  have  a  great  heal  ofelTect  on  the  arms  transfer 
re-u.t  otner  ’han  making  the  growth  of  Evaluate  Negotiate 
reath  c:  the  state  variable  Reject  Refuse  by  about  5  days, 
at  tl.e  diii'erence  m  model  behavior  from  altering  event 


r-.:.;me:ers  s;  ..  '.v  -hat  we  ^annot  maike  any  solid  policy  analysis  recommendations 


tion  aicne. 


B.  CONCLLSIONS 

fl'.e  KSIM  model  has  demonstrated  that  it  is  not  only  feasible,  but  desircable  to 
.ombme  sott  data,  svmh  as  expert  opinion,  with  hard  data,  such  as  the  Third  Point 
S;.  stems  esent  data  set.  'Ehis  capability  to  augment  hard  data  with  subjective  estimates 
s.r;t:.al  '.vl'.en  modelling  systems,  like  the  arms  transfer  process,  which  are  not  well 
vles^rt-'ed  bv  determ.inistic  relationships  with  easily  quantified  parameters. 

.\not;'ier  valuable  implication  of  using  KSIM  is  that  it  provides  a  method  for  the 
ana!;,  st  to  identil'y  the  structure  of  the  system  he  is  modelling,  even  though  he  may 
i'.a'.  e  no  experience  in  modelling  techniques.  The  process  of  estimating  cross  impact 
'..dues  and  event  parameters  will  give  the  user  a  great  deal  of  insight  into  the  "inner 
workings  of  his  sys-.emi.  Often  this  insight  is  just  as  valuable  to  the  analyst  as  actual 
runs  of  the  si.mulation  itself. 

Tile  question  cf  model  validation  is  bound  to  arise  when  discussing  forecasting 
models  such  as  the  KSIM  simulation  of  arms  transfers.  The  answer  is  simply  that  the 
liole  purpose  ci'the  the  model  is  to  provide  the  foreign  policy  analyst  with  a  means  to 
develop  a  mental  image  of  the  arms  transfer  process.  It  is  intended  to  be  neither 
predictive  nor  prescriptivei  merely  a  mathematically  sound  method  for  showing  the 
geometric  relationships  between  system  variables. 

^^T.lle  the  KSI.M  procedure  appears  to  be  a  useful  tool  in  developing  a  model  of 
tile  arms  transfer  process,  the  utility  of  this  tool  in  modelling  various  arms  transfer 
policies  seems  questionable.  The  miodel  is  very  responsive  to  changes  in  the  cross 
impact  parameters,  but  resistant  to  changes  in  the  event  parameters.  However,  the 
various  arms  transfer  policies  are  modelled  by  changing  the  relevant  event  parameters. 
If.us  the  policy  analyst  will  not  be  presented  with  a  clear  picture  of  the  results  of  his 
pvl.c;.  l.ke  those  that  were  seen  when  cross  impact  parameters  were  changed.  Since 
KsIM  relie'.;  on  tiie  geometry  of '.tate  variables,  subtle  changes  in  the  values  of  the 
V  ;:  ;.i'’'le''  cr  time  is  not  sulTcient  to  draw  conclusions  about  the  effect  of  a  given 
.irnt'  'ran-dor  policy. 
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Fieure  4.  /  Results  of  FIxtended  Conaressional  Debate. 
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Some  comment  must  be  made  here  regarding  the  implementation  oi'  a  KSIM 
procedure  with  a  panel  of  experts  with  limited  mathematical  experience,  fhe  ability  to 
understand  the  dilTerence  between  constant  cross  impacts  and  derivative  cross  impact^ 
is  verv  dilTicult  to  someone  who  does  not  possess  a  working  knowledge  of  basic 
calculus.  Thus,  the  policy  analyst  may  easily  becom.e  ballled  when  tr>ing  to  estimate 
parameters  for  these  cross  imipacts  and  when  trxmg  to  understand  why  a  state  variable 
IS  so  aiTected  by  the  rate  of  growth  of  other  variables  .  Furthermore,  the  process  of 
arriving  at  a  final  set  of  cross  impact  parameters  is  one  of  trial  and  error.  This  is  not 
only  very  time  consuming,  but  the  panel  of  experts  have  no  knowledge  of  what  the 
ssstem  behavior  actually  is:  they  only  know  what  they  think  it  should  be  or  uhat  the> 
want  It  to  be.  Clearly,  the  final  arms  transfer  model  could  be  heavily  biased  and  yield 
results  that  are  completely  inaccurate.  To  be  really  useful  as  an  arms  transfer  polic> 
analysis  tool  then,  the  KSI.VI  procedure  needs  some  method  for  obtaining  the  values  of 
the  cross  impact  terms  from  a  hard  data  base. 
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APPENDIX 

KSIM  ARMS  TRANSFER  SIML  LATION  PROGRAM  LISTING 


Programed  in  Turbo  Pascal  (version  3.01a) 


FRCC-RAM  KSIM; 


ParartiRecord  =  RECORD 
IHigh  :  REAL : 

Tlow  :  REAL; 

Max? rob  =  REAL; 

El  ID  ; 

Done , SubHenul , SubMenuZ , SubMenuS :  BOOLEAN ; 
Event , Raw , I , J , RndInitialT , T :  INTEGER ; 
StartTime  ,  StopIi.Tie  .Mini :  INTEGER  ; 

Bij Value , AijValue , Time , DeltaT :  REAL ; 
TiiighVal ,  TlovVal  .MaxProbVal :  REAL  ; 

LcvT ,HighT , MaxP :REAL ; 

Alpha , Beta , InitialT :REAL; 

Name:  SIRING [ 18 ’; 

Bi-Array:  ARRAYf 1 . . 3 , 1 • . 3]  of  REAL; 


£.1  ur  ; 

PlotPoints:  ARRAY [1. .1000,1. ,8]  of  REAL; 
EventData;  AR.RAY  1 1 .  . 3 . 1 .  . 3J  of  REAL; 
iDOintOl:  ARRAY  [l.-Sl  of  REAL; 
EventNames:  ARRAY1I..8]  of  STRING[18]; 
Errorirap:  ARRAYli..8]  of  BOOLEAN; 
Underflow:  ARRAY[l..8]  of  BOOLEAN; 
Cverflov;:  ARR.V/U--8]  of  BOOLEAN; 

Params:  ParamRecord; 

DataOut:  FILE  of  ParamRecord; 

RawData:  FILE  of  ParamRecord; 

Blj File ,Ai: File :  TEXT; 

Names:  TEXT; 

XvsTIME:  TEXT; 

FILEN:  TEXT; 

FUNCTION  PCWER(Term,Exponent:REAL) :  REAL; 

E  C  G  .L  .'I 

POWER  :=  EXP(£xponent*LN(Term) ) ; 

END  ; 

CONST 

Limit:  INTEGER  =  8; 

Empty:  INTEGER  =  -lOCO; 

PROCEDURE  LoadEventNamesArray ; 

BEGIN 

ASSIGN  (Na.mes,  '  Events  .Nam ')  ; 

RESET  (Names;; 

FCR  Event  :=  1  TO  Limit  DO 
BEGIN 

READLN ( Names , Name ) ; 
EventNa.mes  [Event]  :=  Name; 

END  ; 

CLOSE  (Names); 


This  function  computes  the  value 
of  "term"  raised  to  the  power  of 


Loads  an  array  of  the  event 
names  for  use  in  promting 
inputs  and  displaying  some 
of*  the  results  to  the  user, 
Events. Nam  is  on  the  disk. 


%* 


ROCEDURE  LoadData; 

This  procedure  takes  the  estimates  for  times  of  highest  probability 
ar.d  lowest  probabilitv  ^ defined  as  that  point  where  CDF  =  0.10)  and 
the  maximum  probability.  These  values  are  stored  on  the  disk  in  the 
file  RawData.Dat. 

e:-;m 

CIRSCR; 

GOTCX"/  (15,4); 

’;;•■RITEL^;_^  _This  procedure  will  ask  you  to  estimate  parameters'); 

■.-.■RITZi:;  acout  the  events  that  go  into  an  arms  transfer.'); 

R  a  w  D  a  t  a  ,  RawData.Dat'); 

REWROTE  (RawData); 

For  E'.’e.nt  ?=  1  to  Li.mit  DO 
EEDOM 

G0T0XY(24,3) ; 

WRITE  ('For  '  ,EventMames[Event] ,  '  please  enter:'); 

GOTCXrMS  ,  12)  ; 

write  ('The  time  when  the  event  is  most  likely  to  occur  '); 
READLM  (ThighVal); 

Rarams. Thigh  .-=  ThighVal; 

Ever.tData[tvent .  1]  ;=  ThighVal; 

GOTCaY; 15 , 14) ; 

write  ('The  time  before  which  it  probably  will  not  occur  ') 
READLK  (TlowVal;; 

Farams.Tlow  :=  TlowVal; 

Ever.tData[Event,2]  ;=  TlowVal; 

G0T0aY(15,1&) ; 

WRITE  ('The  maximum  probability  of  occurence  '); 

.RE.ADLM  (MaxProbVal)  ; 

Params .MaxProb  :=  MaxProbVal; 

Ever.tData [Event ,  3]  :=  MaxProbVal; 

WRITE  (RawData,  Params); 

CLRSCR; 

EMD  ; 

CLOSE  (RawData); 

W.RITELM(  ‘  Returning  to  Main  Menu'); 

DELAY (500) ; 

eo;d  ; 


ROCEDURE  LoadSioMatrix ; 

This  procedure  takes  the  user's  estimates  of  cross  impact  values.  The  Bi^ 
matrix  is  structured  such  that  IMPACTING  events  are  the  columns  and  IMPACTED 
events  are  the  rows.  Thus  B[3,4]  is  the  impact  of  event  4  cn  event  3. 

The  ^matrix  is  loaded  into  an  array,  BijArray,  and  written  to  the  disk. 

EGO  M 

CLRSCR; 

GOTOXYdO.Z)  ; 

WRlTcLM(  '  This  procedure  will  load  the  B  cross  impact  matrix.  For  each  of) 
GOTOXYdO  ,3)  ; 

WRlTELN('the  following  event  pairs,  estimate  the  cross  impact  of  the  the'); 
GOTOXYdO,4)  ; 

WRITELN( ' first  event  on  the  second.  Positive  values  are  enhancing,  and'); 
GCT0XY(10,5); 

WRITELN( ' negative  values  are  inhibiting.  Use  the  following  scale:'); 
GOTOXY(33,7j ;WRITELN( '0  ...  NONE'); 

G0T0XY(33,8  ;WRITELN('l  ...  MILD'); 

GOTOXY  33,9);WRITELN('2  ...  STRONG'); 

GOT0XY(33,10) ;WRITELN('3  ...  OVERWHELMING'); 

WRITELN; 

ASSIGN  (Bi] File , ' Bi jData . Dat ' ) ; 

REWRITE  di]File); 

FOR  Event  : =  1  to  Limit  DO 
EcGIN 

FOR  Row  :=  1  TO  Limit  DO 
BEGIN 

WRITE  ('The  impact  of  ' .EventNames [Event] , '  on  ’, 
EventNames[Row] , '  ' ) ; 


READi:;  '2i:Value'); 

3i;Array  ]Rov; ,  Event:|_  :=  Bi^Value,- 
See  description  ot  procedure  for  indexino 
values 

V'RITE  I.  3ij  rile  ,  B  i;  Value  )  ; 

WRITEL'.* ; 


EMD  ; 

CIRSCR; 

GOTl'XY^  25 , 2  )  ;WRITELM(  ’0  ...  MOKE  '  1  ; 
jOTCaY  .  25.3;  •'.vRUEiM^^  '1  ...  :!IL2  '  )  ; 

GCTOaYi  25 ,4  1  rWR'TELMi  '  2  ...  STRlMG^r 
G;pAYj_25 , 5  ;  ;V;R;7Elt;(  '  3  ...  OVERWKELMIMG'  )  ; 
V<RI  TELIi ; 


Writes  the  scale 
again  so  it  is 
a -'•■‘a vs  available 


R 1 T  E  L 


(  Gross  i.-ipact  matrix  is  loaded ')  .-WRITELM' ; 
1  Returninc:  to  mam  menu'  '/  ; 


CEGGRE  loadAi^.Matrix : 

■IS  procedure  takes  the  user's  estimates  of  constant  cross  impact  values. 

,e  Ai:  .matrix  is  structured  such  that  IMPACTING  events  are  the  coluiTins 
lu.-.ns  and  I.MF.AGTED  events  are  the  rows.  Thus  A[3,4]  is  the  impact  of  event 
on  event  3 .  .he  matrix  is  loaded  into  an  array,  AijArray,  and  v;ritten 
t  n  e  c  1  s  rw . 


CLR3GR: 
gg::aY"1G.2) : 

wriTaL'.'.' '  Tnis  orocedure  will  load  the  A  cross  impact  matrix.  For  each  of) 

G  vv  1  G  ■ 


V.RITEl.'i;  '  the  following  event  pairs,  estimate  the  cross  i.mpact  of  the  the '  )  ; 
GGTGaYi i; . 4 ; ; 


WRITEEM: • first  event  on  the  second.  Positive  values  are  enhancing,  and'); 
g;t:;<y(ic.5)  ; 


tive  values  are 

inhibiting.  Use 

WRITELN ('0.0  ... 

.  NONE ' ) ; 

WRITELN'/ 'C.l  ... 

.  MILD' ) ; 

■writel.n(  '  0 . 2  ... 

.  STRONG' ) ; 

;WRITELN(  G.3  .. 

.  .  OVERWHELMING'  ) 

ASSIGN  (Ai;File  ,  '  AijData  .Dat '  )  ,- 
REV-'RITE  'Ai:File)r 
FOR  Event  :=  1  to  li.mit  DO 
BEGIN 

FOR  Row  :=  1  TO  Limit  DO 
BEGIN 

WRITE  (‘The  impact  of  ' ,EventNames [Event] , '  on 
EventNames [Row] , '  ' ) ;  ' 

READLN  ^AijValue); 

Ai;Array [Row , Event]  :=  AijValue; 

See  description  or  procedure  for  indexino  of  these 


See  oescription  of  proc( 
values 

WRI.E  .AijFile , AijValue ) 
WRITELN: 


END  ; 

CLRSCR; 

GOTOXY(25,2) ;WRITELNf '0 
G0T0XY(  25,3;  :WRITEi:;(  '  i 
GCT0KY( 25 , 4 }  jWRITELN ' ' 2 
G0T0XY(25,5) ;WRITELN; '3 
WRITELN ; 


NONE '  ) 
MILD'  ) 
STRONG 


OVERWHELMING' 


Writes  the  scale 
aoain  so  it  is 
alwavs  available 


ITEL:.’  (  Cross  imcact  matrix  is  loaded ')  .-WRITELN  ; 
vRITEL;:  -''Returning  to  m.ain  m.enu  ) 


PROCEDURE  Corr.puteKSIMinputs ; 

This  procedure,  reads  tne  raw  data  from  RawData.Dat  and  ccmoutes  Alpha,  Beta 
and  TpcintOl  (the  time  when  the  event  CDF  is  at  .01).  The’  results  are  in 
array  AlphaBeta  and  TpointOl.  AlphaBeta  also  holds  the  value  of  MaxProb. 

Log9:  REAL  =  2.19722457; 

BEGIN 

;-li.nT  :=  ICGO;  Arbitrarily  large  value  so  that  the  minimum  TpointCl  can 
Event  :=  0;  be  found. 

CLR5CR ; 


:o  Limit  DO 

High!  :=  Eventdata [Event , 1  ]  ; 
Lo'.-.’T  :=  Eventdata  [Event ,  2  j  ; 
.MaxP  :=  Eventdata  [Event  ,3]; 


E.ND  ; 

CLOSE  (RawData); 


AlphaBeta [Event , 1 1  :=  Alpha; 
AlphaBeta 'Event, 2'  :=  Beta; 

AlphaBeta [Event , 3]  :=  Max? ; 
TpointOl [Event]  :=  InitialT; 


Inputs  from  RawData.Dat 

:omputation  of  parameters 

)  / ( -Alpha ) : 

)UKD(:nitialT) ; 


Ls  pro 
in  t h ^ 
LiArra 


i  RetrieveBitMatriX; 

Dcedure  recalls  the  cross  impact  matrix  from  the  disk 
;  system  is  ir.tialized  and  builds  the  cross  impact  array 
ay)  for  use  in  the  simulation. 

BIGM  (Bi]File , ' BiiData.Dat ' ) ; 

SET  (BiiFile); 

FOR  J  ;=  1  to  Limit  BO 
BEGIN 

FOR  I  :=  1  to  Limit  DO 
BEGIN 

RE.AD(Bi]File,BijValue) ; 

BijArray [ I , JJ  :=  BijValue; 

See  description  of  the  procedure  LoadBijMatrix 
for  the  indexing  of  this  arrav. 

END  ; 

END  ; 

:sE  (si:?iie); 


SURE  RetrieveAiiMatrix ; 

prccedure  recalls  the  cross  impact  matrix  from  the  disk 
the  system  is  intialized  and  builds  the  cross  impact  array 
irray)  for  use  in  the  simulation. 

ASSIGN  (AijFile, ' AijData.Dat' ) ; 

RESET  (Aljrile); 

FOR  J  ;=  1  to  Limit  DO 
BEGIN 

FOR  I  :=  1  to  Limit  DO 
BEGIN 

READ(AijFile,AijValue) ; 


Ai;]Array[I ,  Jj  :=  AijValue; 

See  description  of  the  procedure  LoadAijMatrix 
for  the  indexing  of  this  array. 


END  ; 
(Ai^File; ; 


PROCEDURE  RetrieveEventData ; 


h: 


1 1.)  in  (.j  1 1  ui:-;  ini  'u  i ;  ;»iui  ''ii.)incii7ii)  <)  ’o  i ;  ri  ■n;iia"  r»  rfu 


rocedure  recalls  the  rav  data  from  the  disk  when  the  crocram 
laliced.  K5I/!  inputs  -alpha,  beta,  and  initial  time)’ are 
ed  i.T, mediately  by  procedure  CcmputeKIShinputs . 


55I3t;  '  R auO a t  a  ,  '  Rawlata  .Cat '  )  ; 

ESET  '  RavC  ata  ' 

IR  Event  :=  1  to  Limit  CC 
EEGi;; 

RE.RC  (  RawData  ,  Params )  ; 
WITH  Params  CC 

EventCata  ^Event , 1 ] 
EventCata [uvent , 2 1 
uventCata[r.vent,  l] 

E^:c  ; 

LCSE  iRawData); 


Params . THioh ; 
Params .  TLov; ; 
Params .  .ha.xProb  ; 


RE  CisplavEventData r 

rccedure  displays  the  current  values  of  event  data  on  the  screen. 


:  Ih’TEGER: 
eturn:  CH.RR  ; 


LRSCR; 

jTOHLWES  ,2;  ;WRITELN( 'E  VENT 
:  T  C  ;-:V  (  2  5  ,  -4 ;  WR I  IE  L!'J  ('Most 
:T0:-;Y(24,5  )  ;WRITELM(  'Likely 
TTCXY  (2S  ,o)  WRITELM  ( '  Time 
IR  Event  ;=  1  to  Li.mit  CO 


PARAMETERS'); 
Not’ )  ; 

Before  Maximum'); 

Time  Probability') 


Column 

headers 


BEC-IM 

V  :=  2’'Event+C; 

G0T0XY('2  ,Y)  ;WRITELN(EventNames  [Event]  ) ;  Row  headers 
GCTC;<:y(25,YK-NRITELN(EventData[Event,lJ  :3:0,  ' 

EventData[Ever.t,2]  :3 :0  '  Curren 

EventData[Event,3]  -A-A)  ;  values 

EMC  ; 

CTC;-:Y^20,24)  ;WRITE('Hit  "ENTER"  to  return  to  main  menu'); 


RE  C isplavBi] Matrix ; 

rccedure  displays  the  current  cross  impact  matrix  on  the  screen. 

1  :  INTEGER  =  6 ; 

2  ;  INTEGER  =  16 ; 


,Y:  INTEGER; 
eturn:  CH.AR; 


;WRITELN('C  ROSS 
;WRITELN( ' IMPACTING 
;WRITELM('l  2 
to  Limit  DO 


I  M  P  A 
EVENTS ' ) ; 

3  4 


BEGIN 

Y  :=  2*I+C1; 

IF  1=4  THEN 
E~GIN 

G0T0XY(5,Y) ;WRITELN( 

END  ; 

IF  1=5  THEM 


T  MATRIX'); 
Column  headers 
5  6  7 


IMPACTED  '  )  ;  ROW¬ 


S'); 


header 


BEGIN 

GOTOXY(6,Y) ;WRITELN( 'EVENTS' ) ;  Row  header 

END; 

GGTCXY ( 1 3 , Y ) ; WRI TELN ( I )  ; 


63 


FOR  J  :=  1  to  Limit  DO 
BEGIN 

X  .=  6*J+C2; 

GOTO:"i(X,Yl ;WRITELN(Bi]Array[I , J] ;2;1);  Current 
EliD;  values 

END; 

GOTOXY(23,24) ;WRITE( 'Hit  "ENTER"  to  return  to  main  menu'); 

READLN  ^Return); 

0; 

OCEDL’RE  DisplavAi jMatrix ; 

'.;is  procedure  Sisplays  the  current  cross  impact  matrix  on  the  screen. 

.'.'57 

Cl:  INTEGER  =  6; 

C2:  IN’TEGER  =  1£  • 

X,Y:  INTEGER; 

Return:  CHAR; 


CLRSCR; 
G0TCaY{23,2; 
ncTryv! Tfi  a 


WRITELN('C  ROSS  IMPACT  MATRIX'); 
WRITELN( ' IMPACTING  EVENTS ' ) ;  Column  headers 
WR:tELN('1  2  3  4  5  6  7  8'); 


c-ctc;-:y(23,5)  ;wr:teln('i  2  3  4  5  6  7  8'); 

FOR  I  :=  1  to  Limit  DO 
BEGIN 

Y  :=  2*I+C1; 

IF  1=4  THEN 
BEGIN 

G0T0XY( 5, Y) ;WRITELN( ' IMPACTED' ) ;  Row  header 

END  ; 

IF  1=5  THEN 
BEGIN 

GOTOXY (6, Y) ;WRITELN(' EVENTS' ) ;  Row  header 

END; 

GOTOXY ( 1 8 , Y ) ; WRITELN (I ) ; 

FOR  J  :=  1  to  Limit  DO 
BEGIN 

X  :=  6*J+C2; 

GOTOXY(X,Y) ;WRITELN(AijArray[I, J] !2:1) ;  Current 
END;  values 

FMO . 

GCT0Xy"23;24) ;WRITE( 'Hit  "ENTER"  to  return  to  main  menu'); 

READLN  (Return) ; 


PROCEDURE  EditEventData; 

This  procedure  allows  the  user  to  change  current  values  of  the  event  data. 
The  current  values  are  displayed  and  new  values  are  input  by  the  user.  The 
new  value  then  replaces  the  old  value  in  the  event  data  array. 

CON’ST 

C:  INTEGER  =  6; 

VAR 

Y:  INTEGER; 

NewValue:  REAL; 

BEGIN 

CLRSCR ' 

GOTOXY(21,3) ;WRITELN('E  VENT  DATA  EDITOR'); 

FOR  Event  :=  1  TO  LIMIT  DO 
BEGIN 

Y  :=  2*Event+C; 

G0T0XY(23,Y) ;WRITELN(Event , ' .  ' ,EventNames [Event] ) ; 

END  ; 

G0T0;"/(  15  ,  (Y+'3 ))  ;WRITE(  '  Enter  the  number  of  the  event  you  want  changed.  ' 
RE.ADLNl Event)  ; 

CLRSCR; 

GCTCXY(22, 3) ;WRITELN( 'Current  values  for  ' ,Eventnames [Event] )  ; 

GCTCXYv 22 , 11 ) ;WRITELN( ' 1 .  Most  likely  time  =  ' , EventData [Event , 1  ]  : 3  : 0 ) ; 


3:T3;<7(22 , 13  );V,'r:tEL.';(  '  2 .  Mot  Before  time  =  ',EventDat3[Event,2;:3:': 

C-CTIXT ;  22 , 15  )  rWRITELla  '  3  .  MaxLrr.'am  probability  =  '  ,  Ever.tData  [ Ever.t ,  3  ^  :  4  : 4 
GCT0XY(  22  19 )  ;V.'R:tE(  '  Enter  the  rinmber  you  vant  to  change. 

READLN  ; I ^  ^  ^  ' 

jOTOXY(  22 , 21 )  r  V.'RITE  (  '  E.nter  the  new  value.  ')• 

RE.ADLH  (MewValue.): 

E'/entData  [Event ,  I  ■  :=  l.’ewValue; 

33T0XY(  22 , 23  ^  ;WR;tEL:.'(  '  Chance  complete  ...  returning  to  .main  menu'); 
CEl.AYdSOO  'i  ; 


;YRE  £dit3i;Matr ix ; 

procedure  allows  the  user  to  change  current  values  of  the  cross 
:t  matrix.  First  the  i.mpacting  event  is  selected,  then  the  i.mpacted 
:  IS  selected.  Che  new  value  is  i.nout  and  repaces  the  old  value  in 
:il.':atrix  array. 


Y:  IMTEGEK; 
h'ewValue:  RE.AL  ; 


31RSCR; 

CCrcXY) 15 , 3) ; 

■.■.■R:rEi:;i  cross  impact  matrix  editor'); 

FCR  Event  ;=1  to  Limit  DO 
5EGIM 

Y  ;=  2*Event+C; 

GOTOXY'v  3C  ,y.)  ;WRITELN(Event ,  '  .  '  .EventNames  [Event]  )  ; 

HMD  : 

GCrc.\Y(2j , ^Y*3) ) ;WRITE( 'Enter  the  number  of  the  IMPACTING  event. 
RE.ADLM  (J); 

CLRSCR- 

FOR  Event  :=1  to  Limit  DO 
BEGIN 

Y  :=  2’'Event+C; 

GCTOXY(30, Y)  ;WRITELN(Event,  '  .  '  ,EventNames[Eventl); 

GOTCXY  20,  )■>  ;WRITE(  '  Enter  the  number  of  the  IMPACTED  event. 

READLN  I); 

CLRS  DR ; 

GCDCXY:  20  , 13 )  ;WRITE  (  Enter  the  new  value  for  the  impact  of); 

GOroXY:  20 , 1 5  ;  .-WRITE  (Event.Names  [  J]  ,  '  on  '  ,  EventN’ames  [  I  ]  ,  '  ' 

READLN  >NewVaiue); 

3i: Arrav[ I , Jl  :=  NewValue; 

GCrc.XY' 2j  .  1=)  ;WRITELN(  '  Change  complete  ...  returning  to  main  menu'); 


:elal  ;1oGu;; 


lEDL'RE  EditAiiMatrix; 

.s  procedure  allows  the  user  to  change  current  values  of  the  cross 


3act  matrix.  First  the  impacting  event  is  selected,  then  the  impacted 
nt  IS  selected.  The  new  value  is  input  and  repaces  the  old  value  in 
Bi] Matrix  array. 


3:  INTEGER  =  6; 

Y:  INTEGER; 

.NewValue:  REAL 

CLRSCR; 

g:t:xy( is  .2) : 

WRIIELi;',  CROSS  IMPACT  MATRIX  EDITOR'); 
FCR  Event  :=1  to  Limit  DO 
BEGIN 

Y  :=  2’^Event+C; 

G0TCXY(30,y) ; WRITELM(Event , ' .  ' , EventNames [Event] ) ; 


END  • 

GOTOXY  V  20  (Y-^3) )  ;WRITE(  '  Enter  the  number  of  the  IMPACTING  event. 
READLN  (J): 

CIRSCR; 

FOR  Event  :=1  to  Limit  DO 
BEGIII 

Y  :=  2^Eventi-C; 

GOTC:<Y V  3 ,  Y )  ;WR:TELN(Event,  '  .  '  ,  Eventh'ames  [ Event  ]  )  ; 

END  ; 

GOTO;<Y :  2  0  .( Y+3  ));  WRITE  (' Enter  the  number  of  the  IMPACTED  event. 


GDTCaY  !  20 , 1 3 )  rWRITE  (  '  Enter  the^new  value  for  the  im.pact  cf‘); 
GCTOXY ;  20 , 1 5 )  ■  WRITE  ^Ever.tNames  [  J]  ,  '  on  '  ,  EventMam.es  [  1 1  ,  ' 
RZADLM  (NevValue^- 


Ai  jArrav^I,  Jl  :=  MevValue,- 

GCTOXYt  ul ,  13)  ;WRITEL.‘I(  ■  Change  comolete  ...  returning  to 

w  -  -.A  -  \  ^  D  ^  ' 


main  menu 


PRCviEYRE  InitializePlctPointsArrav ; 

Tne  arrav  PlctPoints  holds  the  values  of  each  event  at  each  time  increment. 
BEGIN 

FOR  T  :=  0  to  1000  do 
BEGIN 

FOR  Event  .-=  1  to  Limit  DO 

PlotPoints [T .Event]  :=  0; 

END  ; 

IT'”'  , 

PROCEDURE  InitializeTime ; 

This  procedure  fi.nds  the  start  time  for  the  simulation.  It  selects  the 
minimum  TpointOl  unless  that  value  is  less  than  0.  If  this  occurs  then 
the  start  time  is  set  at  0  and  CDF  values  are  computed  for  those  events 
with  Tcoir.tOl  less  than  0. 


Ir.itializePlotPcintsArrav : 
FCR  Event  ;=  1  to  Lim.it  ic 


BEGIN 


ErrorTrap [Event]  :=  FALSE; 
Underflcvi Event j  :=  FALSE; 
Overflov[Ever.t]  :=  FALSE; 


Xvalue [Event , ll 
.Xvalue  [Event  ,2] 


=  0;  Initializes  Xvalue  array.  Xvalue [t,l] 
=  0;  is  X(t  *  dt).  Xvalue[t,2]  is  Xit). 


CL.RSCR  • 

GGTOX'/ f  17  , 12 )  .-WRITE (  ' .How  many  days  do  you  want  to  simulate?  f250  max) 

READLN  StcpTi.me )  ; 

DeltaT  :=  0.25; 

FOR  Event  -.=  1  to  Limit  DO 
BEGIN 

Alpha  :=  AlohaBeta [Event, 1] ; 

Beta  :=  AlphaBeta [Event , 2 j ; 

MaxP  :=  AlohaBeta^Event , 3] ; 

IF  TpointCl [Event]  <  0  THEN 
BEGIN 

ToointOl [Event]  :=  Empty; 

Xvalue [Event , 1 '  :=  MaxP/ ( 1  +  EXP(Beta));  CDF  at  t=0 
Xvalue [Event , 2 ]  :=  MaxP/(l  +  EX?(Eeta));  CDF  at  t=0 
Co.m.outes  X(0)  since  T(.Ol)  <  0 

END  ; 

END  ; 

IF  Mini  <  0  THEN  StartTime  :=  0  ELSE  StartTime  :=  HinT ; 

G jT0XY( 30 , 16 ) ;WRITELN) ^ Start  time  will  be  '.StartTime); 


PROCEDURE  I terativeEqua tion ; 


\  a> 


*  J"  r  V,  V.  tr-.  « 


.nis  procedure  computes  the  value  of  X( t  +  dt)  based  on  the  value  of  /hti. 

It  t3.<es  into  account  cross  irr.pacts.  The  equations  are  described  in  CHATTER 
3  or  the  thesis. 

AR 

GarrjnaSunPlus  , Gao'”.aSu.r.Minus  :  REAL; 

Gam.v.a  .  PIsubI  ,  dHdt .  Isno.T'. , h’urr.  :  REAL; 

X.plusDeltaT  :  REAL: 

Aloha  :=  AlphaBeta[I  1’; 

Seia  :=  AlphaE^-a'  l  2’ '■ 

Gav-aSuoPlus  ;=  0- 
Gar'-aSunHinus  :=  C  ■ 


SEGi;; 


IF  ! /lvalue  [  J  ,  i :  =  Xvalue[J,2j)  THE.'T  Gam.’na  :  =0 


BEGIM 


dXdt  :=  'Hivaluel. 


1  ;  -  .'.va^ue 


^/.  t  ^  +  d  u  \  /  /  C  t 

Gair..~ia  :  =  .Ai::.Ar  rav[  I  ,  J  j  ^Xvalue  1  J  ,  1  ’  -  Bi' Ar  rav  1 1  ,  J  '  ''dXd 

EMD  : 

Ganv.aSuriMinus  :=  Gatnr.aSu.T.Minus  +  ABSlGam.va)  -  Gam.va; 

Ga.T-.aSu.v.Flus  ;=  GaT,ma5u.T.r  lus  t  AB5  (Garri.T,a ;  -i-  Gao..r.a; 

END ; 

Pisubl  :=  (1  +  C .  5’'DeltaT’'Gamn'^a5u.TiHinus  )  /  ( 1  +  0 . 5''DeltaT’'Ga~.i',aS'u.'nPlus  )  ; 

IF  Error. rap[Ij  THEH  prevenos  underflov  error  vioh  EXP  ni  Ham  and  Denom 
EE  Gill 

Hum  : =  0 . C  ; 

3-1*0 p.  *“  i.C* 

E1ID_ 

BEGIM 

Mur,  :=  Alpha*Deltat’^EXP( -Alpha’^Time  +  Beta); 

Denom  :=  1  +  EXP( -Alpha''Time  +  Beta); 

IF  {EXP( -AlDha’'Time  +  Beta)  <  lE-20)  THEN  Underflow  error 
Errorirapfl]  :=  TRUE;  is  im.renent. 

END ; 

IF  Underflow!' I]  THEM  XTolusDeltaT  ;=  0.0 
ELSE 

IF  Overflow! I]  THEM  XTplusDeltaT  :=  1.0 

5  “  Q  T  *.J 

XTclusDe:tal‘'7=“pcWER(XvalueLl,ll  ,  PI  subi ) ’'POWER)  ( 1  +  Mum, /Denom)  ,PIs’ubl!  ; 

IF'  IXvaluefl,!]  <  lE-20)  THEM  Underflowill  :=  TRUE: 

IF  'Xvalue[l,lj  >  0.999999)  THEN  Overflow) I]  :=  TRUE ; 

EMD  ; 

Xvalue)I,2]  :=  Xvalue[I  1);  Replaces  X(t)  with  old  X(t  +  dt) 

Xvalueil.ll  ;=  XTolusDeltaT;  New  X(t  +  at)  for  ne.xt  ti.me  increment 
PlotPoints  [T  ,  I  j  Xvalue[I,2]; 

PROCEDURE  Simulation; 

The  simulation  sets  X(t)  =0  if  the  time  is  before  T(.Ol)  for  X.  When  time 
qual  to  or  greater  to  T(.OOl)  then  Xft)  is  set  to  0.001.  Thereafter,  the 
edure  IterativeEquation  computes  values  for  Xit  +  dt) . 

Zero:  RE.AL  =  0.0; 


Ti.me  :=  StartTim.e; 

T  :=  1  ; 

'WHILE  Ti.me  <=  StopTixe  DO 

G:TCX’/(35  ,  IS  )  ;WRITELM(  '  Time  =  ■,Time;3:l); 
.-or  I;  =  l  to  Limit  DO 


IF  {  Tooint  31 '  I  j ''=Time}  .AMD  ( TpointOl  [  I  ]  oEmpty )  THEM 
Checks  if  T).u1;  is  less  than  present  tim.e.  If  this  value 
IS  less  than  present  time  OR  if  it  is  empty,  these  steps  are 

S  hi  ICC-i  d  . 


BEGIN 


few-Jw' 


F? 


I 

A: 


If  »(".  ■. 


•ft? 

m 

L-i: 


B.' 


en: 


Xvalueri,li 
Xvaluefl  2l 
Tpcir.cOl  [Ij 


0.001  ; 
0.001; 
Empty ; 


:ase 


For  I:  =  l  to  Li.mit  DO 
BEGIN 

IF  y.value[I,l]  =  Zero  THEN 
BEGIN 

I  of 

1  :  Plotpoi.''.tS 
2:  Plotpoir.ts 
3:  Plotpoints 
4;  Plotpoir.ts 
5:  Piotpoi.its 
d:  PlotD01."tS 
7  :  Plotpoi.nts 


Both 
i.n  I 
Sets 
will 


END; 


8:  Plotpoints 


are 
terat 
EMPT 
be  sk 


0 

0; 


set  at 
iveEoua 
Y  value 
ipped . 


C.OCl  s 
tio.i  wi 

so  above  c.iec)< 


0  switch 


0 

C 

0 

0; 

0; 

0; 


END 

ELSE  IterativeEquation ; 


END  : 


EN 


END; 

Time  :=  Time  +  Delta! ; 
T  :=  T  +  1 


Increments  time 

Increments  PlotPoints  array  row  index 


u  ; 


P.POCED’JRE  PrintResults  ; 

Sends  results  to  the  printer 
BEGIN 

CLR5CR; 

G0TCKi’(25 , 15 )  ;WRITELN( '  Sending  results  to  printer.') 
DEL.RYf^OOh- 

WRITE ^LST, 'Time  1234 

WRITELN(LST, ' 5  6  7  8'); 

Ti.me  ;=  StartTime; 

WHILE  Time  <=  StopTime  DO 


.  ime 
BEGIN 


[T.l 

:  4  :  4  ,  '  '); 

'T,2 

:  4 : 4 ,  '  '); 

'T,3i 

:  4  : 4  ,  '  '  )  ; 

■T,4j 

!  4  : 4  ,  '  '); 

’T.S 

:  4  :  4  ,  '  '  )  ; 

'T,6 

:4:4,'  '); 

[T;Z] 

:4:4,'  '); 

tr 


BEGIN 


WRITE  (LST,Time:4:l, ' 

WRITE  (LST.PlotPoints 
WRITE  (LST, Plotpoints 
WRITE  (LST.PlotPoints 
WRITE  (LST.PlotPoints 
WRITE  (LST.PlotPoints 
WRITE  (LST.PlotPoints 
WRITE  (LST.PlotPoints 
WRITELN  (LST.PlotPointstTiS] :4:4) 
Time  :=  Time  +  8''DeltaT; 

T  :=  T  +  8; 

END; 

GOTOXY(30.17) ; 

WRITELN( ' Enter  To  Continue'); 

READLN ; 


CEDURE  WriteOutputFile ; 

is  procedure  writes  a  sequential  data  file  to  the  disk.  The  user  inputs 
e  name  of  the  file.  The  procedure  writes  time,  XI  values,  X2  values,  etc. 
me  is  multiplied  by  100  and  rounded,  X  values  are  multiplied  by  10000  and 
unded.  This  is  so  that  the  data  is  not  in  exponential  notation.  The  data 
n  then  be  plotted  on  high  resolution  plotters  such  as  GRAFSTAT. 

Column , RndPoint , RndTime :  INTEGER ; 

Point:  REAL; 

GutFileMame:  STRING[i4]  ; 


CLP.SCR; 
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GOTOrf  (25 , 14)  ;WRITE  (.'What  name  do  you  vant  for  the  output  file' 
GCTOX’/<,2  5,.17); 

READLM  ^OutFiieMame ' ; 

ASSIGN  'XvsTIME,  CutFileN’am.e )  ; 

REWRITE  (MvsTIME); 

Tim.e  :=  Startlim.e; 

REPEAT 

RndTime  :=  ROUND ( Tim.e'^lCG )  ; 

W'RITELM  ( y.vsTINE  ,  RndTirne  >  ; 

Ti.TiS  Ti.t.s  Dsl-sT* 


Time  :=  StartTime; 

repeat' 

Point  :=  PlotPoints [T , Column] ; 
RndPcint  :=  ROUND{Point*10000) ; 
■WRITELM  (KvsTIKE,  RndPoint); 
Time  ;=  Time  +  DeitaT; 

V  T  :=  T  -  1  r 

UNTIL  Time  >  StcoTime; 

EN’D  : 

CLOSE  ;XvsTlNE ) ; 


RCCEDURE  DisplavNainNenu ; 


SucNenul  ;=  FALSE;  Preset  to  FALSE.  Selection  of  a  menu  choice 
Suc;-!enu2  :=  FALSE;  1,2  or  3  will  switch  and  select  proper  s'ubmenu 
Sub.NenuS  :=  FALSE;  interpretation. 

CLRS  CR ; 


WRITELN  ('CROSS  -  INPACT  ANALYSIS  FOR  ARMS  TRANSFERS') 
WRITELN  ( 'HAIM  HENU ' ) ; 

;WRITELN  ('1.  Enter  New  Data'); 

•  wr’'T":m  ('0.  Review  Current  Data' 

Edit  Current  Data ' ) ; 

Run  Model ' ) ; 

Write  Output  File  for  GRAFSTAT'); 

Exit  Program '  ) 


•RITE..N 
;WRITELN 

;Writel;; 

:  'fi'RIZELN 
;WRITELN 
.•WRITE  (' 


)  ; 


o . 

(■  '4. 


T^-pe  n'umber  would  you  like,  then  hit  ENTER 


RCCEDURE  D 1 sp 1 a y S ubN e nu ; 


CLRSCR; 

GCTONY  32,a):WRITELN  ('SUB 
GOTCXY(25 , 10) ;WRITELN  ('1. 
GOTCXY(25 , 12  , -WRITELN  ( ' 2. 
GCTCXY(25 , 14) ;WRITELM  ( ' 3. 
GOTCXY(25 , 16) ;WRITELN  ( '4. 
G0TCXY(17 ,21)  , -WRITE  ('T’/De 


MENU ' )  ; 

Event  Parameters'); 

Constant  Cross  Impact  Matrix'); 
Derivative  Cross  Impact  Matrix'); 
Return  to  Main  Menu '  )  ,- 
number  would  you  like,  then  hit  ENTER 


RCCEDURE  ExitProgram; 

This  program  writes  the  current  values  of  the  cross  impact  matrix  and 
the  event  data  to  the  disk  and  exits  the  program. 


CLRSCR; 

ASSIGN  (BijFile  ,  '  BijData .  Dat '  )  ,- 
REWRITE  (BijFile); 

FOR  J  :=  1  TO  Limit  DO 
BEGIN 

FOR  I  :=  1  TO  Limit  DO 
BEGIN 

BijValue  :=  Bi jArray [ I , J] ; 
WRITE  (BijFile ,Bi jArray [I , Jj)  ; 


END  ; 


END ; 


CLOSE  (BijFile)r 

ASSIGN  (AijFile, ' AijData .Dat ' ) ; 

REWRITE  \AljFile); 

FOR  J  :=  1  TO  limit  DO 
BEGIN 

FCR  I  :=  1  TO  Limit  DO 
BEGIN 


AijValue  :=  AijArray[I , J] ; 
WRITE  ( Ai] F ile , Aij Array [ I , J] 


END: 

END  ; 

CLISE  (AiiFiie); 

ass:  on  Rai'Data  ,  '  RawData  .Dat '  )  ; 

REWRITE  ;RavData): 

FOR  I  ^=^l^TO  Limit  DO 

Rarams.THigh  :=  EventData [ I , 1 ] ; 
Params.TLov  -.=  £ventData[1 , 2]  ; 
Para.T.s  . MaxProb  :=  EventData[1 ,3]  ; 
WRITE  ( RawData , Params ) ; 

E;:D; 

^-^OSr.  (RawData)  : 

:-;to:-:7;so  .  i4)  ;Writeln(  -  end  of  analysis'  ) ; 

DELAY i 1500) ; 

CLRSOR; 

Dtr.e  :=  TRYE  ; 


lURE  ErrorFlag; 

procedure  is  activated  if  the  user  selects  an  invalid  choice  from 
r.er.u.  It  always  returns  to  the  main  menu. 

WRITELN  ('INVALID  MENU  SELECTION  ...  CHOOSE  AGAIN'); 

DELAY  (7C0);  Allows  user  to  read  message. 


SURE  InterpretMainOption; 
rprets  menu  selection  from  the  main  menu 

Option  :  INTEGER; 

Response  :  CH.AR; 

READLN  (Option); 

CASE  Ootion  OF 
1  :  BEGIN 

SubMenul  :=  TRUE;  Insures  proper  sub  menu  actions 
DisplaySubMenu; 

END  ; 

2:  BEGIN 

SubMenu2  :=  TRUE; 

DisplaySubMenu; 

END  ; 

3;  BEGIN 

SubMenu3  :=  TRUE; 

DisplaySubMenu ; 

END  ; 

4:  BEGIN 

ComputeKSIMinputs :  Computes  event  parameters 
InitializeTime ; 

Simulation;  Actual  KSIM  model  in  here 

G0T0XY(23, 20) ;WRIT£( 'Print  Results?  y  for  YES 
READLN( Response ) ; 

IF  Response  =  'y'  THEN  PrintResults ; 

END; 

5:  WriteOutputFile ; 

6:  ExitProgram; 

ELSE  ErrorFlag;  Returns  to  the  main  menu 


-EDURE  InterpretOptionl ; 

Ls  procedure  activated  vhen  #1  is  the  main  menu  choice.  These 
rform  the  sub  menu  actions  requested. 

SubOption  :  IMTEGER; 

::i 

RE.iTLM  (SubCption)  ; 
r.i.SE  SubCDtibn  OF 
1  :  ‘  EEC-IM 

LoadData ; 

Com.DuteKSIMinputs  ; 

EMD  ; 

2:  LoadAiijMatrix; 

3:  LoadBi:Matrix; 

4:  5EGIM 

CLRSCR; 


END; 

E13E  ErrorFlag; 


GOTGXY(  30 , 15 )  ;V;RITELN  (' Returning  to  Main  Menu’) 
DELAY  (500); 


Returns  to  the  main  menu 


•RCDED'JRE  InterpretOptionD  ; 

This  procedure  activated  when  #2  is  the  main  menu  choice.  These 
_cispj.ay  the  data  requested. 

SubOption  :  INTEGER; 


RE.nDLN  ^SubOption) 
C.RSE  SubCotion  CF 


DisplayEventData 
Display.-.ij  Matrix 
DisplavBijMatrix 
BEGIN 


CLRSCR • 

GOTOXY( 30 , 15) ;WRITELN( ' Returning  to  Main  Menu') 
DELAY  (800); 


ELSE  ErrorFlag;  Returns  to  the  main  menu 


lED'JRE  Interpret0ption3  ; 

13  procedure  activated  when  #3  is  the  main  menu  choice, 
rform  editing  on  the  data  base  requested. 

SubCption  :  INTEGER; 


These 


?::.iDLN  (SubOption); 

SASE  SubOption  OF 

1:  ‘  EditEventData ; 

2:  EditAiiMatrix; 

3:  EditBi]MatriX; 

4:  BEGIN 

CLRSCR; 

GOTCXY(30, 15) ;WRITELN( 'Returning  to  Main  Menu') 
DELAY  (600); 

END; 

ELSE  ErrorFlag;  Returns  to  the  mam  menu 
EL'j  ; 


■:a::;  frcgram 


TEXTBACKGROUND  Q) 
TEXTCCLCR  (14;,- 


Done  :=  False; 

RetrieveEventData ;  Initializes  the 

RetrieveBijMatr IX ;  system  with  the 

RetrieveAijMatrix;  current  values. 

LoadEventNamesArray ; 

WHILE  NOT  Done  Do 
BEGIN 

DisplayMainNenu; 

InterpretHainOption; 

IF  SubMenul  ..HEN  InterpretOptionl  ; 
IF  Submenu2  THEN  Interpretoption2 ; 
I?  Submenus  THEN  InterpretCptionS ; 

END  ; 
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